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SUMMARY 


Under contract NASl-13655, Grumman developed the computer code CAVE 
(Conduction Analysis via Eigenvalues) , a convenient and efficient computer code for 
predicting two-dimensional temperature histories within thermal protection systems 
for hypersonic vehicles. NASA report CR-2897 describes fully the CAVE code and its 
operation. 

Under the present contract the capabilities of CAVE have been enhanced by 
incorporation of the following features into the code: 

• Real-gas effects in the aerodynamic heating predictions 

• Geometry and aerodynamic heating package for analyses of cone-shaped 
bodies 

• Input option to change from laminar to turbulent heating predictions on 
leading edges 

• Modification to account for reduction in adiabatic wall temperature with 
increase in leading edge sweep 

• Geometry package for two-dimensional scramjet engine sidewall, with an 
option for heat transfer to external and internal surfaces 

• Print -out modification to provide tables of select temperatures for plotting 
and storage 

• Modification to the radiation calculation procedure to eliminate temperature 
oscillations induced by high heating rates. 

This report describes these new features and is an addendum to report CR-2897. 



Section 1 


INTRODUCTION 

The CAVE code (Reference 1) , which was developed for LaRC under contract 
NASl-13655, predicts the transient temperature response within two-dimensional 
geometries. It has found application in transient heating studies of thermal pro- 
tection systems for hypersonic aircraft and cruise missiles. CAVE was designed 
with the convenience of the user in mind. Usual operation of CAVE requires only 
such minimal information from the user as material properties, initial temperature 
distribution, flight trajectory (i. e. , altitude, velocity, and angle -of-attack as 
functions of time in tabular form) and selection of one of the built-in general struc- 
tural configurations, with overall dimensions and grid network sizes specified via 
input data. CAVE discretizes the geometry into elements and calculates the 
capacitances, conduction links, and convective boundary conditions. The transient 
temperature response of the structure is calculated using a hybrid analytical- 
numerical technique in which spatial derivatives are replaced by appropriate finite 
difference representations and the temporal derivatives are retained as ordinary 
derivatives. The numerical technique, developed under LaRc contract NAS-1-11818 
(Reference 2) , is inherently stable and permits large time steps through the flight 
trajectory with attendant savings in computer time. 

Under the present contract, the utility and capabilities of CAVE have been 
enhanced by incorporation of the following features into the code; 

• Real -gas effects, which become important for Mach numbers above six, into 
the aerodynamic heat transfer predictions 

• A geometry and aerodynamic heating package for analyses of cone-shaped 
bodies (Figure 1) 

• Geometry package for two-dimensional scramjet engine sidewall, with an 
option for heat transfer to external and internal surfaces (Figure 2) 

• Input option to change from laminar to turbulent heating predictions on 
leading edges 
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• Modification to account for reduction in adiabatic wall temperature with 
increase in leading edge sweep 

• Print -out modification to provide tables of select temperatures for plotting 
and storage 

• Modification to the radiation calculation procedure to eliminate temperature 
oscillations induced by high heating rates. 

As part of the present contract, a Thermal Analyzer Network Generator was 
provided to LaRc. This code will be used in conjunction with CAVES (Conduction 
Analysis via Eigenvalues for 3-dimensional Geometries) which was developed under 
contract NASl-14643. CAVES (Reference 3) is a general computer code for analysis 
of problems involving any geometry, boundary conditions, and materials. This 
generality is achieved by requiring the user to establish the thermal network, e.g. , 
node capacitances, conductances between nodes, etc. The Network Generator will 
significantly reduce the manpower expended creating and debugging models. Basically 
the Generator utilizes a digitizing tablet for geometry data input, model generation 
software for calculation of such parameters as node capacitances and conductances, 
and graphics hardware for displaying the model. A separate report discusses the 
details and operation of the Network Generator. 

The present report is an addendum to the original user’s manual for CAVE 
(Reference 1). Section 2 discusses the modifications that have been incorporated to 
account for real -gas effects which become important at Mach numbers above 3 or 4. 
Section 3 discusses the new Blunt-Nosed Cone Geometry; and Section 4 discusses the 
new Scramjet Engine Side-Wall Geometry. For both geometries, a sample problem is 
given in addition to the required input data format. 

In Section 5, the new option for creating special output files of temperatures and 
heat transfer coefficients at pre-selected nodes is given. This file can be useful for 
making quick comparisons and for generation of time-history plots. Section 6 describes 
the modification that has been added to analyze Leading -Edge Geometry problems with 
laminar to turbulent flow transition. Section 7 outlines the modifications that have 
been incorporated into the technique for handling radiation heat transfer. These 
modifications are directed at eliminating the temperature oscillations that have on 
occasion been induced by heat heating rates (i. e. , rapid changes in temperature) and 
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large time steps. The Appendix presents flow charts for the new subroutines that 
have undergone more than trivial modification. 

As part of the present contract, the CAVE code was carefully reviewed relative to 
the reduction in adiabatic wall temperature with increase in leading-edge sweep. 

Careful examination of the code indicates that the reduction in T^^ with increase in 
sweep angle is accounted for fully and no changes were required. 

Mr. James L. Hunt of the High Speed Aerodynamics Division, Langley Research 
Center, Virginia, served as the NASA technical monitor for the program. 

At Grumman, the contract was administered by the Advanced Development 
office under Mr. Fred Berger, Manager of Advanced Development System Engineering. 
Dr. Kenneth A. Rathjen was Study Manager and Mr. Henry O. Burk was Programmer. 
The many helpful discussions with Dr. Robert L. Kosson are gratefully acknowledged. 
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Section 2 


REAL-GAS EFFECTS 

For Mach numbers less than 3 or 4, there is virtually no error in calculating 
the total temperature of a gas via the equation: 

Tt = Ts + (1) 

where = total temperature, R 

Tg = static temperature, R 
"t = c /c 

p V 

M = Mach number of gas 

This equation is valid for gases that satisfy the perfect gas law (Pv = RT) , have con- 
stant specific heats, and undergo an adiabatic process. For Mach numbers above 4, 
the temperature dependency of the specific heats introduces errors in the calculation 
of total temperature via Eq,. (1). At Mach 6, the error is approximately 10%, and at 
Mach 7 it is approximately 20%. This error is on the order of 400 to 500 °R, which 
can be quite significant in the evaluation of thermal protection systems for hypersonic 
vehicles. 

To extend the usefulness of CAVE, the method of calculating the freestream 
total temperature was modified to the following for Mach numbers greater than three: 

1. For the given altitude find, via subroutine ATMOS, the atmospheric ambient 
pressure and temperature T^. 

2. For the given velocity and ambient temperature calculate the freestream 
Mach number: 

Mee=Xo/(32.2TRTJ^/^ 

3. Calculate the dynamic pressure: 

2 



4. Compute pitot pressure: 

Pt 2 8* 92 X 10~^ q^i^atm 

-4 

where the constant 8. 92 x 10 is equal to the product of 14. 7 x 144 x K 
with K = 0. 53. 

5. Compute freestream enthalpy: 
h«,= C g J T or 

h^=0.24 • 25036 • ft^/sec^ 

6. Compute a normalized total enthalpy: 

(h/RT^) = (0 . 5V^ + hj /843914 

t ^ 2 

where the constant 843914 is equal to the product of g RT with 

c o 

R = 53. 3 ~ Ip and T =491.6 R= 32 F. 
ibi - R o 

7. Via a table look-up, using (h/RT ) and log (P ) establish the freestream 

O , 2 -to t, ^ 

total temperature T. * 

t, z 

The table used in Step 7 is to be the first table in the data set of any run. It is 
provided as part of the code and is referred to as the Standard Table. The format for 
the table is that required by subroutine NUREDl described in Appendix A of 
Reference 1. Figure 3 displays the 46 card images of this table. 
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2714 


TEMPERATURE AS A FUNC OF LOG 1 01 PRESSURE) S ENTHALPY OOTTGT 
7*6 7*723 B.42 9*12 9«7a5 10*485 11*191 13*343 14*080 OlTTOT 

-2.0 1059*fi5l080« 1170* 1260* 1350* 1440* 1530* 1800* 1890* 02TT0T 

-1 *69901 OS9*85lOeO* 1 170* 1260* 1350* 1440* 1530* 1800* 1890* 03TTOT 

-1*3^791 059* £5 1080* 1 170* 1260* 1350* 1440* 1530* 1800* 1890* 04TTGT 

- 1 * l£491059.eS10eO* 1170* 1260* 1350. 1440* 1530* 1800* 1890* 05TTOT 

-1*0 1059*851080* 1170* 1260* 1350* 1440* 1530* 1800* 1890* 06TT0T 

-0*69901059*851080* 1170* 1260* 1350* 1440* 1530* 1800* 1890* 07TTOT 

**0«397910S9.CS1G80» 1170* 1260* 1350* 1440* 1530* 1800* 1890* 08TTOT 

-0*1 £491059*851080* 1170* 1260* 1350* 1440» 1530* 1800* 1890* 09TTQT 

0*0 1059*851080* 1170* 1260* 1350* 1440* 1530* 1800* 1890* lOTTGT 

0«30 1031059.651080. 1170. 1260. 1350* 1440. 1530. 1800* 1890* IITTOT 

0*602061059*851080* 1170* 1260* 1350* 1440* 1530* 1 799* 301 889*951 2TTOT 

0*84£101059*£51080* 1170* 1260* 1350* 1440* 1529* 1799* 1 889*4913TTGT 

1.0 1059.851060* 1170. 1260* 1350* 1440* 1 528.75 179 6. 76 1689* 1 01 4TTOT 

1*30 1031 059*651080* 1170* 1259* 1350* 1433*581527*551798* 151 887*601 5TTOT 

16*351 19*362 22*520 25*801 29*4 33*905 37*39 40*2 42*5 1 6TTGT 

-2*0 2160* 2520* 2880*553240*0 3600* 3960* 4185* 4310* 4400* 17TT0T 

-1*69902160* 2520* 288 1 • 123242* 093609*4540 03* 7842 1 0* 4365* 4470* 1 8TTOT 

-1*39792160* 2520* 288 1 • 193242*6 136 1 6* 0 740 18* 644260* 4440* 4550* 19TTGT 

-1*1 £4921 60* 2520* 286 1 *553244 • 233620 *324036* 854300 * 4495* 4600* 20TTOT 

-1*0 2160* 2520* 2881*803245*2 3622*684047* 144320* 4505* 334637* 452 1 TTGT 

-0*6 9902160* 2520* 2881 * 803245* 753625*5 1 4062* 674350*75455 1 • 234695*6022TT0T 

-0*39792160* 2520* 288 1 • 803246* 253627 *8 74076* 084379 *474595*674754 .S523TT0T 

-0* 15492160* 2520* 288 I • 803246* 273629*2940 85* 194394*524625* 084793* 7224TTGT 

0*0 2160* 2520* 268 1 *803246* 2 73629*7840 88* 544407*71 4640* 54461 5* 192STTOT 

0*30 1032i59*952S19*952eai*803246.273630*7|4094*924422*a94665* 744£4S*9526TTOT 
0*60 2062 159* £5251 8* 83268 1 *353246* 2 7363 1 * 1 84099* 804 434*924687* 384680* 4927TTGT 
0* 8451 021 59 *6025 18*24286 1* 133246* 25363 1 * 1 84 1 0 1 * 624 441 * 75470 1 * 27490 1 *2 126TTQT 
1*0 21 58*86251 7*952880*663245* 7 53630*88 4 I 02 • €54 445 *994707* 7449 1 0 • 7729TT0T 

1*30 1032 157* 7225 I 7*662879*093244* 253629 *884 103* 284448*9147 17* 074926* 1 130TTOT 
45*0 50*0 60*0 70*0 60*0 90*0 100*0 110*0 120*0 31TTOT 

-2*0 4490* 4670* 4945* 5170* 5435* 5790* 6325* 6905* 7550* 32TT0T 

-1*69904585* 4765* 5065* 5310* 5575* 5925* 6400* 7120* 7665* 33TTOT 

-1*39794675* 4870* 5175* £465* 5745* 6065* 6530* 7195* 999994-934TTOT 

-1*15494720* 4945* 5290* 5570* 5665* 6185* 6610* 7260* 999994 - 93STTOT 

-1*0 4760* 4985* 5340* 5645* 5930* 6260* 6690* 7305* 999994-936TTGT 

-0*69904820* 5060* 5455* 5600* 6110* 6485* 6£45* 7435* 999994937TTQT 

-0*39794915* 5185* 5595* 5945* 6265* 6615* 7025* 7560* 999994-938TTOT 

- 0 * 15494975* 5230* 5700* 6065* 6410* 6745* 7 160* 7680* 999994-939TTGT 

0*0 5000* 5285* 5740* 6145* 6515* 6880* 7295* 999994-9999994940TTOT 

0*30 1035035* 5360* 5885* 6310* 6675* 7050* 7480* 9999949999994-94 ITTOT 

0*602065070* 5420* 5995* 6460* 6850* 7255* 7580* 99999 49999994-942TT0T 

0.84510S095 * 5450* 60 9 5* 65 8 0 * 7010 * 74 25* 9999 94-9999994-9999994-943TTQT 

1*0 5175* 5505* 614E* 6630* 7085* 7515* 999994-9999994-9999994-944TTOT 

1*30 103515 0 * 5520* 6200* 677 0 * 7260* 999994-9999994-9999994-9999994-945TTOT 


1827-003B 

FIG. 3 TABLE OF 2 AS A FUNCTION OF (h/RTQ>^ 2 & LOG (P^ 2> 


9 




Section 3 


BLUNT-NOSED CONE GEOMETRY 


3.1 DISCUSSION 

This subsection presents the blunt-nosed cone geometry that has been 
incorporated into CAVE and discusses how this geometry is discritized into nodes 
by the BLUNT2 subroutine. Subsection 3.2 presents the input data format for this 
geometry and Subsection 3. 3 presents a sample problem of a hollow cone. 

Figure 4 shows a cross-sectional view of the blunt-nosed cone geometry assumed 
by subroutine BLUNT2 when it generates a nodal network. The insulating layer can 
be eliminated by using an input value of zero for TAU. A solid cone is selected by 
using JGEO = 4, and a hollow cone is selected by using JGEO = 5. Optionally, the 
hollow cone may be cooled internally. Figure 5 shows the grid network that is gener- 
ated by BLUNT 2 for solid and hollow cones. Because CAVE is a two-dimensional 
code, only cones at zero angle -of-attack can be analyzed. Consequently, it is suffi- 
cient to show only the upper half of the cone with the understanding that the nodal 
elements shown are to be rotated about the axis of the cone to produce a body of 
revolution. 

Unlike the cooled leading edge geometry, the cooled cone geometry places no 
restrictions on the choices for AX^. The thickness of the insulating material is given 
by TAU and it may equal zero. Notice in Figure 5 that nodes are located on the inter- 
face between the insulating material and the main cone material. BLUNT2 assumes 
that there are equal volumes of the two materials associated with each interface node, 
i. e. , that one half of ^^ 3 * this case, is associated with the insulator and the other 
half with the main material. 

The user may elect to have CAVE calculate the convective heat transfer 
coefficients and adiabatic wall temperatures over the cone or he may supply tabular 
inputs for them. If the user selects the former option, tabular values for the flight 
parameters of velocity and altitude (angle -of-attack must equal zero) as functions of 
time must be supplied. The code uses the axisymmetric method of Fay and Riddell 
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(Reference 4) for the stagnation point heating value and Lee’s method (Reference 5) for 
the heat transfer distribution on the spherical nose. Downstream of the spherical nose 
the Mangier transformation factor is applied to obtain the proper heat transfer 
coefficients for the cone. 

Problems involving increased heating due to local interference heating or 
other effects such as plume impingement during part of the trajectory can be handled 
by using the HMODI flags in the same way they are applied in the leading-edge 
geometry case. It is noted that the surface node elements modified in this way 
correspond to an entire "ring" around the cone and not a small local surface area on 
one side of the cone. 

In using this multiplier option, it is important to bear in mind that in modifying 
the convective coefficients, CAVE takes the average of the multiplicative factor at 
the beginning and end of the time interval and applies it over the entire interval. 
Therefore, the tables and computing times intervals must be selected with some care 
whenever step changes are to be simulated. A sample problem in Section 4 of 
Reference 1 illustrates this point. 

3.2 INPUT DATA FORMAT FOR BLUNT-NOSED CONE GEOMETRY 
Table Output Card 

• ITPRNT (II) 

ITPRNT = 0: Standard data tables will not be printed 
ITPRNT = 1: Standard data tables will be printed at beginnii^ of 
printed output. 

Standard Tables Cards 

The Standard Table is read in so that the program can perform table look-ups 
of temperature versus enthalpy and pitot pressure. This method is used to obtain 
TTOT at Mach numbers >3.0 (thereby accounting for real gas effects). The table 
must always be read in as part of the input data package whether it will be used or not. 

The input format for this table is described in the section dealing with Sub- 
routine NURED, given in Appendix D, Sheet D-1 of Reference 1. 
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• Table 1 contains: (46 cards) 

TTOT (R) as a function of ) 

° t,2 

(Arg. 1) and Log^^^ (P^ g) 2) , where ^ i® *^^e pitot pressure 

(in atmospheres) 

• Blank Card (signifies end of standard tables). 

This table is provided as part of the computer deck. 

Indexes Card 

• JGEO, L, M, NE (415) 

JGEO = 4 (selects solid cone geometry) 

= 5 (selects hollow cone geometry) 

L = number of elements through the material (may be 

an even integer or odd integer) 

M = number of elements along top surface of cone 

NE = number of dominant eigenvalues to be used in solution 

(e.g. , a typical number is 5). 

Title Card 

• Run identification, comments, etc. (8A10) 

Radiation Card 

• EPSl, TBGl (2F10.5) 

EPSl = emissivity of surface 

TBGl = background radiation temperature, R. 

Material Properties Cards 


• MAT (15) 

• NMATl, RHOl, CONAVl, CPAVl (I10,3F10.5) 

• TMATl(l), TMAT1(2), . . . , TMATl (NMATl) ^ omit (8E10.0) 

t CONDTl(l), CONDTl(2) CONDTl (NAMTI) > if (8E10.0) 

• CPTl(l), CPT1(2) CPTl (NMATl) 1 NMAT1=0 (8E10.0) 


(If MAT = 2 include the cards:) 



• NMAT2, RH02, CONAV2, CPAV2 

• TMAT2(1), TMAT2(2) TMAT2 (NMAT2) 

• CONDT2(l), CONDT2(2) CONDT2 (NMAT2) 

• CPT2(1), CPT2(2) CPT2 (NMAT2) 


(I10.0.3F10.5) 
omit (8E10.0) 

if (8E10.0) 

NMAT2=0 (8E10.0) 


MAT 

NMATl 

RHOl 

CONAVl 

CPAVl 

TMATl(l) 

CONDTl (I) 
CPTl (I) 


= number of materials (1 or 2) 

= number of entries in properties table (maximum 
of 10) . NMATl = 0 for constant properties 
= density of material 1, Ibm/cu-ft 
= average thermal conductivity of material 1 (used 
when NMATl = 0) , Btu/ft-sec-°R 
= average specific heat of material 1 (used when 
NMATl = 0), Btu/lbm-°R 
= temperatures in thermal properties table for 
which CONDTl (I) and CPTl (I) are given; 

1=1, 2 NMATl, °R 

= thermal conductivity of material 1 at temperature 
TMATl (]), Btu/ft-sec-°R 
= specific heat of material 1 at temperature 
TMATl (I), Btu/lbm-°R 


NMAT2, RH02, CONAV2, etc., same as NMATl, RHOl, CONAVl, etc., 
except applied to material 2 


Detail Geometry Cards 

• MCAP, THETA (110, F 10. 5) 

• DELX(l), DELX(2), DELX(3) DELX(MM) (8F10.5) 

• DELR(l), DELR(2), DELR(3) DELR(L) (8F10.5) 

• TAU (Flo. 5) 

• RP, HCOOL, TCOOL (3F10.5) 

(omit this card when GEO = 4) 


MCAP = number of elements into which nose of cone 

is subdivided (must be an even integer) 

THETA = wedge half angle of cone, in large 

DELX(I) = spatial increments in X direction 1=1, 2 MM 

(where MM = M - MCAP) , ft 
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TAU 


= thickness of material 1, ft (when considering only one 

material, TAU = 0) 

RP = radius of nose coolant passage, ft 

HCOOL = convective heat transfer coefficient inside one coolant 

2 

passage, Btu/ft -sec- R 
TCOOL = cone coolant temperature, R 


Initial Temperature Cards 

• KOBE, I, T(I), II, JJ (215, ElO. 0,215) 

• ... 

• ... 


11100 (indicates end of initial temperature cards) 

KOBE = 0 or blank 

I = node number 

T(I) = node initial temperature, R 

II and JJ = the node number is incremented by the spacir^ JJ until 

the limit II is reached. Each node so specified is 
assigned the same temperature. 


(15) 


Special Output File Card(s) 


• NPLOTS, ITCOBE (214) 

• NOBNUM(l), NOBNUM(2), ..., NO BNUM (NPLOTS) 

(omit if NPLOTS = 0) 

NPLOTS = total number of nodes to be placed in special output file 

ITCOBE = 1 if node temperatures are to be stored in R 

= 2 if node temperatures are to be stored in K 
NOBNUM(I) = node number of temperature to be stored 1=1, 2, ... 
NPLOTS. 


Boundary Condition Cards 

Two options exist: 1) the user inputs the flight trajectory and the code calculates 
the convective boundary conditions along the top surface of the cone and 2) the user 
inputs directly the convective heat transfer coefficient and adiabatic temperature as 
functions of time and distance. 
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OPTION 1. FLIGHT TRAJECTORY SPECIFIED 


REFX, CODEX, HMODI 
GAM, RGAS, PR 
NT RAJ 



(3E10.0) 


(3E10.0) 


(110) 

TIMTAB (NTRAJ) 

(8E10.0) 

ALTTAB (NTRAJ) 

(8E10.0) 

VELTAB (NTRAJ) 

(8E10.0) 


REFX = not applicable to cone geometry 

CODEX = 0 (indicates to code that Option 1 is being exercised) 

HMODI = nonzero value indicates that a table will be read at 

the end and used to multiply the convective coefficients 
GAM = ratio of specific heats of air 

RGAS = gas constant for air, ft-lbf/lbm-R 

PR = Prandtl number of air 

NTRAJ = number of points in trajectory table (2< NTRAJ <50) 

TIMTAB(I) = time in trajectory table 1=1, NTRAJ, sec 

ALTTAB(I) = altitude corresponding to time TIMTAB(l) , ft 

VELTAB(I) = velocity corresponding to time TIMTAB(l) , ft/sec. 

OPTION 2. CONVECTIVE COEFFICIENT AND ADIABATIC WALL 
TEMPERATURE SPECIFIED 




REFX, CODEX, 
REFX 
CODEX 
HMODI 


HMODI (3E10.0) 

0. (not applicable to cone geometry) 

-1. (indicates to code that Option 2 is being exercised) 
nonzero value indicates that a table will be read at the 
end and used to multiply the convective coefficient. 


Convective Coefficient & Adiabatic Wall Temperature Tables 

Two tables are required. The first gives the convective coefficient as a function 
of time (argument 1) and distance (argument 2) . The second table gives the adiabatic 
wall temperature as a function of time (argument 1) and distance (argument 2) . In 
setting up these tables. The tables must be followed by a blank card. The specifics 
on the format for the tables are given in the description of subroutine NUREDl in 
Appendix A of Reference 1. 
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The following input data is required for both options: 

Time Intervals Cards 

• NTIMES (110) 

• TIMES(l), TIMES(2), ... TIMES (NTIMES) (8E10.0) 

NTIMES(I) = number of points in time intervals array (2 < NTIMES < 50) 

TIMES(l) = initial time (usually equals 0.) sec 

TIMES(I) = time at which temperature is calculated and printed 
1=2, 3 NTIMES, sec. 

Convective Coefficient Modification Table (omit when HMODI = 0.) 

One table is required to modify the convective coefficient on the aerodynamic ally 
heated surface. The multiplicative factor is given as a function of time (argument 1) 
and distance (argument 2) . The writeup for subroutine NUREDl (Appendix A of Ref- 
erence 1) gives the specifics on the format requirements. Follow this table with a 
blank card; omit the blank card if the table is not read in. 

3.3 SAMPLE PROBLEM FOR CONE GEOMETRY 

This subsection contains an illustration of the solid blunt-nosed cone geometry 
considered as a sample problem. Figure 6 shows the grid network. The main fea- 
tures are that it is made of aluminum, has a nose radius of 0.4 ft, and has a length 
of 1. 75 ft. The boimdary conditions of h and uniform and constant with time. 

The following pages show listings of the input and output data for this problem. 
The sequence of the output is: 

• Statement regarding storage allocated for S array in main program 

• Comment card 

• Volumes of nodes computed by CAVE 

• Geometry related input data 

• Node numbers adjacent to exterior boundary 

• material properties 

• Boundary conditions table 
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• Node number location within output array 

• Material number assigned to each node. (In this problem there is only one 
material being used.) 

• The capacitance of each node 

• The conductance in the x-direction between nodes 

• The conductance in the y-direction between nodes 

• Initial temperature distribution. 

And then the following information is printed for each time interval: 

• Average heat transfer coefficients calculated using the h values at the 
beginning and end of this time interval 

• Average heat transfer couplings which include radiation effects , if any 

• Average adiabatic wall temperatures for this time interval 

• Temperatures at the end of the interval 

• Steady-state temperatures for the heat transfer couplings and adiabatic 
wall temperatures of this interval 

• Integrated heat input to each node. This gives the net heat transfer at 
each boundary node up to the end of this time interval. 

(Annotation has been added to the input and output to aid the reader.) 
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NOTE: ITPRNT, STANDARD TABLE & BLANK CARD NOT SHOWN. 


J680 L 
1 1 

M NE 
i i 


T T 

4 S 

T T 

15 5 


TEST CASE FCR 

SOLID BLUNT 

NCSED 

0.0 0.0 



0 

169. O.0AJ06 

s 

10. 


0.05 

0. 1 

0.2 

0.2 

0.2 


o.os 

0.1 

0.1 

0.0 

^ ’ TAU 


0 1 

SOD. 75 

‘1 



-I— 


COMMENT CARD 
NO RADIATION 
1 MATERIAL 
AVERAGE MATERIAL PROPERTIES 
MCAP, THETA 

0.2 0.2 0.2 0.2 


■AR'S 
VTURE 

SPECIAL OUTPUT FILE OF SELECT NODES, IN R 
REFX, CODEX, HMODI 


-AX'S 


02 06 


blunt 6CSE0 CCNE HT OIST. 

00 


0.0 

155. 

01 

0.0 

.00555 

•00555 

02 

0.3 

.00555 

•00555 

03 

.75 

•005ES 

•00555 

04 

1.50 

•00555 

•00555 

OS 

2.25 

•00555 

•00555 

06 

2.7 

•00555 

•00555 

07 

0206 


BLUNT NOSED CCNE TAH DIST. 

00 


0.0 

155. 

01 

0.0 

1500. 

1500. 

02 

0.3 

1500. 

1500. 

03 

• 75 

1500. 

1500. 

04 

1. SO 

1500. 

1500* 

OS 

2.25 

1500. 

1500. 

06 

2.7 

1500. 

1500. _ _PI pAOn 

07 


TABLE OF H VALUES 
FOR 2 TIME VALUES 
& 8 DISTANCE VALUES 

TABLE OF Taw VALUES 
FOR 2 TIME VALUES 
& 8 DISTANCE VALUES 


0. — TIME STEP INTERVALS 


NOTE: THE TIME 
& DISTANCE VALUES 
HAVE TO BE THE SAME 
IN THESE TWO TABLES 
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SHEET 3.3 INPUT/OUTPUT FOR BLUNT-NOSED SOLID CONE PROBLEM 
bD (REFER TO FIG. 4) (SHEET 1 OF 6) 

CO 




4' t- t 4 • 


a' 






15 SURFACE ELEMENTS — 


V 

y 

5 ROWS BT 


• CAVE CODE 

* 

* 4 44444444 44 444444444 * 

IS CCLUMNS GIVES 
















0 








75 ELEMENTS 


y' 

5 DOMINANT MOOES. • .REQUIRES 


ECCNCNOMIZE...REOUCE DIMENSION OF S AND VALUE OF MhOROS FROM 11000 TCHAROS 


2301 MQR05 OF MEMORY 

VALUE REQUIRED FOR THIS PROBLEM 


4 444444444 4 «444*4444444* 444*4444444*44444 44* 444*4 444 444*444*4 44* 444* 444*4 44* 444*444*444*444*44444444444«444«4444444A4444 

TEST CASE FOR SOLID BLUNT NOSED CONE 


VOLUMES OF NODES 
















V( 

11= J.0017139 

VI 

2)> 

0.0022109 

VC 

31 = 

0.0009939 

VC 

*) = 

0.0002637 

VC 

5) = 

0.0000101 

VC 

6) = 

0.0050090 t 

V( 

7)= 0.006*613 

VC 

8)» 

0.00290*6 

VC 

9)» 

0.0007706 

VC 

10) = 

0.0000296 

VC 

11) = 

0.0079160 

VC 

121 = 

0.0102112 1 

VI 

13)= 0.00*590* 

VC 

1*) = 

0.0012178 

VC 

15) = 

0.0000*68 

VC 

16) = 

0.0102097 

VC 

17) = 

0.0131699 

VC 

161 = 

0.0059204 1 

VI 

191= 0.0015707 

VC 

20) » 

0.000060* 

VC 

211 = 

0.01 1712* 

V< 

22) = 

0.01S10B3 

VC 

23) = 

0.0067916 

VC 

2*) = 

0.0018019 I 

VI 

251- 0.0000693 

VC 

261 = 

0.0058732 

VC 

271* 

0.0047130 

VC 

28) s 

0.0031660 

VC 

29) = 

0.0016191 

VC 

30> = 

0.000*569 f VOLUME OF 

VI 

31 )= 0.0122196 

VC 

32) = 

0.0098992 

VC 

33)* 

0.006005* 

VC 

3*) = 

0.003711 5 

VC 

351 = 

0.0013911 

VC 

36) = 


VI 

37 >= 0.0220763 

VI 

38) = 

0.0158806 

VC 

39) = 

0.0097009 

VC 

*0) = 

0.0050601 

VC 

*1 )■ 

0.0300535 

VC 

*2) = 

0.025*127 ) NODES CALCULATED 

VI 

*3>= 0.0192250 

VC 

**) = 

0.0130373 

V( 

45)* 

0.0003965 

VC 

46) = 

0.0339030 

VC 

*7) = 

0.0292622 

VC 

*8) = 

0.02307*5 1 nw oAWrr 

VI 

*9)= 0.0160866 

VC 

50) = 

0.0122*60 

VC 

51 ) = 

0.036265* 

VC 

52) = 

0.03362*6 

VC 

53) = 

0.027*369 

VC 

S*) = 

0.0212*92 1 B* CAVE 

VC 

5S)= 0.016606* 

VC 

56) = 

0.0*31 409 

VC 

57) = 

0.0385001 

VC 

58) = 

0.032312* 

VC 

59) = 

0.02612*7 

VC 

60) = 

0.021*839 1 

VI 

611= 0.0*0529* 

VC 

62) = 

0.0*36006 

VC 

631 = 

0.0377009 

VC 

64) = 

0.0315132 

VC 

65 ) = 

0.026872* 

VC 

66> = 

0.05**31)8 1 

V( 

67)= 0.0*97901 

VI 

68 )> 

0.0436023 

VC 

69) = 

0.037*1*7 

VC 

70) = 

0.0327739 

VC 

71)* 

0.0608*5* 

VC 

721* 

0.05620*6 f 

VI 

731= 0.0500169 

VI 

7*)» 

0.0*38292 

VC 

751* 

0.0391884 

VC 










4 4 44 4444 44 4*444*4 44* 444* 444*444*444*4 44*4 44*444*4 44*444*44444 44*444*444*4444444*444*444*444*444*444* 444*444*444*444*444* 



B L 

U N T NOSE 

: D CCNE 

P R 0 e 



NOSE RADIUS* 

0.4000E400 

FT 


length of 

CCNE SECTICN- 

0.1750E401 

FT 



TMETA* 

0.1000E402 

DEG 



TAU = 

0.0 

FT 



EM1SS1V1TY= 

0.0 



RAOIAT ICN 

6ACKGRCUN0 T= 

0.0 

OEG.R 

the 

EXTERIOR BOUNDARY NODES ARE 




6 

11 16 21 26 31 36 

* 1 *6 51 

56 61 

66 71 


4 44* 444* 4 4 4* 444* 4 4 4* 4 44* 444* 4 4 4*4 4 4* 4 44*4 44* 444* 4 4 4* 4 44«4 4 4*4 44 *444* 44 4* 4 44* 444 *4 44*4 44 *444*444*4 44* 444*444*4 44 *4 44*4 44* 

MATERIAL PROPERTIES CONSTANT 

material 1 RHO=109»00 LBM/CU-FT K=0«OA306 BTU/SEC-FT-OEG.R CP^O.ZOBO BTU/LBM-OEG.R MATERIAL 

PROPERTIES 

4 4 4* 444* 4 44*444*4 44* 4 44* 444*444*4 44*4 44*4 44* 444* 4 4 4*>*44 *444*4 44*444* 444*4 44* 444*444*444*444*444*444«444«444*444*444*444* 
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SHEET 3.3 INPUT/OUTPUT DATA FOR BLUNT-NOSED SOLID CONE PROBLEM 
(REFER TO FIG 4) (SHEET 2 OF 6) 






TABU 

E S 







TABLES 2 AND 3 GIVE H AND TAW 

FOR E)<TERIOR 

BCUNOARV 

NODES CF CCNE 




2 6 

BLUNT NOSED CCNE 

HT OIST. 






TABLES 

0*0 

0*0 0«1590E«>03 0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0*5950E-02 0.9950E>02 0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0.0 

READ IN 

0«3000£«00 

0.9950E-02 0»9950E>02 0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*7S00E «00 

0*5550E-02 0.S5S0E-02 0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

FORh 

O.l 530E»01 

0 .SSSOE- 02 0 • SSSOE- 02 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*2250E«^01 

0*S5S0E-02 0»95S0E-02 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8>Taw 

0.2700E4-01 

0.9S90E-02 0.59SOE-02 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Z 6 

BLUNT NOSED CONE TAW 

OIST. 







0*0 

0.0 0. ISS0E«03 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0*0 

0. 1500E404 0. 1900E404 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0>3C00E4-00 

0.I900E404 0* kSOOE^-OA 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0«7500E»00 

0.150DE404 0.1900E404 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0*1 SOOE^Ql 

0.I900E404 O.1SOOE404 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0»2 2S0E*0I 

0.I900E404 0.1900E404 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0«2700E^01 

0.1900E4-04 0.1S00E»04 0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0 0 









1 



NODAL 

N E T « 0 R K 






hOCE NUMCERS 








ROM / COL 

12 3 

4 

5 

6 

7 

8 

9 

10 

1 

L.O 6.0 IL.O 

L6.il 

21.0 

26.0 

31.0 

36.0 

41.0 

46.0 

2 

2.0 7.0 12.0 

1 7.0 

22.0 

27.0 

32.0 

37.0 

42.0 

47.0 

3 

3.0 B.O 13.0 

18.0 

23.0 

28.0 

33.0 

38.0 

43.0 

48.0 

« 

4.0 9.0 14.0 

19.0 

24.0 

29.0 

34. 0 

39.0 

44.0 

49.0 

5 

5.0 lO.O 15.0 

20.0 

25.0 

30.0 

35.0 

40.0 

45.0 

50.0 

Ron / COL 

11 12 13 

14 

15 






i 

51.0 96.0 61.0 

66.0 

71.0 






2 

52.0 57.0 62.0 

67.0 

72.0 






3 

93.0 Sa.O 63.0 

68.0 

73.0 






« 

94.0 59.0 64.0 

69.0 

74.0 






5 

95.0 60.0 GS.O 

70.0 

7S.0 






NATERI AL 

AUMBER AT EACH NODE 








RQI / COL 

12 3 

4 

S 

6 

7 

8 

9 

' 10 

1 

1.0 1.0 1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

2 

1.0 1.0 1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

3 

1.0 1.0 1.0 

1.0 

1 .0 

1.0 

1 . 0 

1.0 

1.0 

1.0 

4 

1.0 1.0 1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

5 

1.0 1*0 1.0 

1.0 

1.0 

l.O 

l.O 

l.O 

l.O 

1.0 

ROM / COL 

11 12 13 

14 

IS 






1 

1.0 1.0 1.0 

1.0 

1.0 






2 

1.0 1.0 1.0 

1.0 

1.0 




ONE MATERIAL 

3 

4 

1.0 1.0 1.0 

1.0 1.0 1.0 

1.0 

1.0 

1.0 

1.0 




THROUGHOUT 

S 

1.0 1.0 1.0 

1.0 

1.0 






CAPACITANCE AT EACH NODE 








RO« / COL 

1 2 3 

4 

5 

6 

7 

e 

9 

10 

1 

0.«02E>0I 0. 176E400 0.27BE4-00 

0.399E400 

0.4I2E400 

0.206E400 

0.430E4-00 

0.939E+00 

0.106E40I 

0.1 19E401 

2 

0.777E-01 0.227E400 0.359E400 

0.463E400 

0.531E400 

0. 166E-I-00 

0.348E4-00 

0.776Et00 

0.893E400 

0.103E+01 

3 

0.349E-01 0. I02E400 0. 161E400 

0.208E400 

0. 239E4Q0 

0. 11 IE-i-00 

0.239E400 

0.559EA00 

0.676E400 

0.811E400 

4 

0.927E-02 0.27IE'01 0.420E-01 

0.9S2E-01 

0.633E-01 

0.569E-01 

0. 130E400 

0.341E400 

0.4S8E400 

0.S94E+00 

5 

0.356E-03 0.104E-02 0. 16SE-02 

0.212E-02 

0. 244E-02 

0. I61E-01 

0.489E-01 

0. 178E4-00 

0.295E400 

0.430E400 

1827’009B 










SHEET 3.3 INPUT/OUTPUT DATA FOR BLUNT-NOSED SOLID CONE PROBLEM 
(REFER TO FIG. 4) (SHEET 3 OF 6) 




RON / COL 


1 

2 

3 

« 

S 


li 

0* 13SE4-0I 
0.1 18E4-01 
0*964E4>00 
0.747£4^00 
0>Sa4E«>00 


12 

0.IS2E401 

0»135E40t 

0.114E«^01 

0*918E+00 

0.7SSE400 


13 

0« 171E4-01 
0* 1S4E*0I 
0« 133E40I 
0« 1 1 1E401 
0»94SE»00 


14 

0* 19IE401 
0. 17SE401 
0* IS3E401 
0.132E401 
0* 1 15£«-01 


15 

0*214E401 
0» I98E*0I 
0* 176E401 
0* I54E»01 
0* 138E4^01 


conductance in X-DIRECTION. 


□ h 

/ COL 

1 

2 

3 

1 


0.265E>01 

0.190E-OI 

0.307E-0I 

2 


0.S29E'01 

0.379E-01 

0.615E-01 

3 


0.529E-0I 

0.379E-01 

0.515E-01 

4 


0.529Eo01 

0.379E-'0a 

0.61SE-0I 

5 


0.265E-0I 

0.190E-01 

0.307E-0I 

DN 

/ COL 

1 1 

12 

13 

1 


0.944E-01 

O.I01E400 

0. 107EfOO 

2 


0.623E-01 

0.670E-OI 

0.717E-01 

3 


0.489E-01 

0.S36E-01 

0«S83E-01 

4 


0a356E-0i 

0«403E*0I 

0* 4S0E'0I 

S 


0.392E'02 

0.450E^02 

0.509E-02 


0»399E>01 

0.798E-0I 

0*798E"01 

0*798E'01 

0.399 E '-01 

14 

0*1 13E»0Q 
0«764E-01 
0.630E-01 
0«497E>01 
0.568E-02 


0* IOOE+00 
0« 15bE«-00 
0« I26E«00 
O . aOlE -01 
0> 191E-04 
15 

0.0 

0.0 

0.0 ^ — 

0*0 

0.0 


6 

0.686E-01 
0* 1 lOEf^OO 
0.749E-01 
0.394E-01 
0.300E-02 


0.709E-01 

0.587E-01 

0.410E-01 

0.232E-01 

0*209E-02 


e 

0*7S6E-01 

0«482E>01 

0.34BE-01 

0.21SE-01 

0.21SE-02 


9 

0«819E-01 

0.529E-0I 

0.395E-01 

0«262 E ''01 

0a274E-02 


10 

o*aeiE-oi 

0.S76E-01 

0*442E-01 

0*309E»01 

0«333E-02 


ZERO VALUES BECAUSE THERE ARE 
NO NODES TO THE RIGHT OF THESE 


COhOUCTANCE IN VoOlRECTIGN 


RON / CCL 
1 
2 

3 

4 

5 

RON / COL 
1 
2 

3 

4 

5 


1 

0*126E-01 

0.629E-02 

0.210E>02 

0*349£-03 

0.0 

11 

0»26SE4-00 
0.212E400 
0.158E»00 
0. I49E«>00 
0.0 


2 

0*107E4>00 

0.537E-01 

0.179E-01 

0*298E>02 

0*0 

12 

0*284C«00 
0.231E400 
0. 177E400 
Oa 17SE400 
0.0 


3 

0. 170E«>00 
0aB49E-01 
0.283E'01 
0*472E-02 
0*0 

13 

0«303£4-00 
0.2S0E+00 
O* 196E«00 
0.200E4^00 
O.O 


0.219E400 
Oa 109E400 
0.385E-01 
0.606E-02 
Q.Q 

14 

0.322E400 
0a2E9E400 
0.2 1SE400 
0.22SE«^00 
0.0 


S 

0. 251E400 
Oa 126E4-00 
0«419E-01 
0.696E-02 
0«0 


6 

Oa 469E-01 
0.33SE*01 
0* 198E-01 
Oa lOSE-01 


0.973E-01 

0.705E-01 

0.433E-01 

0.263E-01 


8 

0.209E400 
0. 155E + 00 
0*101E»00 
0.725E-01 


9 

0.228E400 
0. 174E4-00 
0. 120E4-00 
0.984E-01 


15 

0.341E4Q0 

0.287E*00 

0.234E+00 

0a2S0E»00_ 

0- 


0 


ZERO VALUES BECAUSE THERE 
ARE NO NODES BELOW THESE 


10 

0.247E+00 

0.193E<»00 

0.139E4-00 

0.124E+00 

OaO 


IN ITIAL 

temperature 

DlSfRIBUTICN 

OEG.R 








RON 

/ COL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 


500.0 

500*0 

500.0 

500.0 

500*0 

500.0 

500*0 

500*0 

500.0 

500.0 

2 


500.0 

500.0 

500.0 

500*0 

500.0 

500*0 

500*0 

500*0 

500*0 

500*0 

3 


500*0 

500*0 

500*0 

500*0 

500*0 

500*0 

500*0 

500.0 

500*0 

500*0 

4 


500.0 

500*0 

500*0 

500*0 

500.0 

500.0 

500.0 

500.0 

SOO.O 

500.0 

5 


500.0 

500*0 

500*0 

500*0 

500*0 

500*0 

500*0 

500.0 

500*0 

500.0 

RON 

/ COL 

1 1 

12 

13 

14 

IS 






1 


500*0 

500*0 

500*0 

500*0 

500.0 






2 


500*0 

500*0 

500*0 

500*0 

500*0 






3 


500*0 

500*0 

500.0 

500*0 

500*0 






4 


500*0 

500*0 

500.0 

500.0 

500.0 






5 


500*0 

500*0 

500.0 

500*0 

500.0 
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4-*4«* T I M E = 0*2000E4-02 SECONDS 4>4>4«« 


AVERAGE 

HEAT TRANSFER 

COEFFICIENTS BTU/SEC-FT**2-0EG.R 







RDM 

/ 

COL 

1 

2 

3 

4 

s 

6 

7 

B 

9 

10 

1 



Q.555E-02 

0.555E-02 

0.555E-02 

0.5S5E-02 

0.55SE-02 

0.SS5E-02 

0.S55E-02 

0.S55E-02 

0«SS5E'02 

0.555E-02 

2 



0*0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

3 



0*0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

4 



0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

5 



0*0 

0.0 

0.0 

o 

• 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROM 

/ 

COL 

11 

12 

13 

14 

15 






1 



0.55SE-02 

0.555E-02 

0.S9SE-02 

0.S5SE-02 

0.S55E-02 






2 



0*0 

0.0 

0.0 

0.0 

0.0 






3 



0.0 

0*0 

0.0 

0.0 

0.0 






4 



0.0 

0*0 

0.0 

0.0 

0.0 






5 



0.0 

0.0 

0.0 

0.0 

0.0 







AVERAGE 

HEAT TRANSFER 

COUPLINGS 

BTU/SEC-DEG.R 







ROM 

/ 

COL 

I 

2 

3 

4 

5 

6 

7 

e 

9 

10 

1 



0.216E-03 

0.632E-03 

0.S9BE-03 

0* 129E-02 

0. I48E-02 

0.694E-03 

0. 143E-02 

0.30SE-02 

0.329E-02 

0.353E-02 

2 



0.0 

0.0 

0.0 

0*0 

0. 0 

0*0 

0.0 

0*0 

0*0 

0.0 

3 



0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0*0 

0.0 

0*0 

0*0 

4 



0.0 

0.0 

0*0 

0*0 

0.0 

0*0 

0.0 

0.0 

0*0 

0*0 

5 



0.0 

0*0 

0*0 

0*0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

ROa 

/ 

COL 

11 

12 

13 

14 

15 






1 



0.37BE-02 

0.402E-02 

0.426E-02 

0.450E-02 

D.475E-02 






2 



0.0 

0.0 

0.0 

0.0 

0.0 






3 



0.0 

0*0 

0*0 

0. 0 

0*0 






4 



0.0 

0.0 

0.0 

0.0 

0.0 






S 



0.0 

0.0 

0.0 

0.0 

0.0 







AVERAGE 

ADIABATIC MALL 

TEMPERATURE 

OEC.R 








ROM 

/ 

COL 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 



1500.0 

1500*0 

1500.0 

1500.0 

1500*0 

1500.0 

1500.0 

1500*0 

ISOO.O 

1500.0 

2 



0.0 

0.0 

0.0 

0.0 

0*0 

0*0 

0*0 

0*0 

0*0 

0*0 

3 



0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0.0 

4 



0*0 

0.0 

0.0 

0.0 

o 

• 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

5 



0.0 

0*0 

0*0 

0.0 

0*0 

0*0 

0*0 

0.0 

0*0 

o 

. 

o 

ROM 

/ 

CGL 

11 

12 

13 

14 

15 






1 



1500*0 

1500*0 

1500*0 

1500.0 

ISOD.O 






2 



0*0 

0*0 

0.0 

0*0 

0*0 






3 



0*0 

0*0 

0.0 

0.0 

0*0 






4 



0*0 

0*0 

0.0 

0.0 

0.0 






S 



0.0 

0*0 

0.0 

0*0 

0*0 
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TEMPERATURES AT THIS 

TINE DEG.R 









ROM 

/ COL 

1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

1 


529* 1 

527*0 

S26.6 

526. 6 

527.6 

529.0 

531.2 

531.2 

530.4 

529.6 

2 


522.3 

522.7 

522.7 

522.6 

523.1 

522.0 

522.7 

S21.9 

520.7 

519.7 

3 


519.0 

519.4 

519.6 

519.5 

519. 1 

516.0 

516.7 

5 14*6 

513.0 

511.7 

A 


517*0 

5 17. 1 

517.1 

516.7 

516* 1 

514.7 

512.5 

509*5 

507. 1 

505.4 

S 


516.3 

516.3 

516.3 

516.2 

516*1 

512.4 

510.4 

506.9 

504,2 

502.5 

ROk 

/ COL 

11 

12 

13 

14 

IS 






1 


S20.a 

527.9 

527.0 

526.4 

526*1 






2 


sia.e 

517.8 

516.0 

516.2 

515.0 






3 


510.5 

509*3 

50 8.3 

507.6 

507.2 






4 


504. 1 

502.0 

501.6 

501.0 

500.7 






5 


S01.2 

499.8 

490.6 

490.0 

497.0 






STEADY-STATE 

TEMPERATURES FOR THE 

BOUNDARY CCNDITICNS AT 

THIS TIME 

DEG.R 





ROM 

/ COL 

1 

2 

3 

4 

5 

6 

7 

e 

9 

10 

1 


1 500. 0 

1500.0 

1 500.0 

1500.0 

1500.0 

1500.0 

1500*0 

1500.0 

1500.0 

1500.0 

2 


1500.0 

1500*0 

1500.0 

1500* 0 

1500*0 

1500*0 

1500*0 

1500.0 

ISOO.O 

1500.0 

3 


1500.0 

ISOO.O 

1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

1500*0 

1500.0 

1500.0 

4 


1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

1500.0 

ISOO.O 

S 


1500*0 

1500*0 

1500.0 

1500.0 

1500*0 

1500*0 

1500*0 

1500.0 

1500.0 

1500.0 

ROW 

/ COL 

1 1 

12 

13 

14 

15 






1 


1500.0 

1500.0 

ISOO.O 

1500.0 

1500*0 






2 


1500.0 

ISOO.O 

1500.0 

1500.0 

1500.0 






3 


1500.0 

ISOO.O 

1500.0 

1500.0 

1500.0 






4 


ISOO.O 

1500.0 

1500.0 

1500.0 

1500.0 






5 


1500.0 

1500*0 

1500*0 

ISOO.O 

1500*0 






INTEGRATED 

HEAT INPUT 

AT EACH NODE BTU 








ROW 

/ COL 

1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

t 


4.2 

12.3 

19.4 

25* 1 

26*7 

13*5 

27.8 

59.1 

63.8 

68.6 

2 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 


0.0 

0.0 

0.0 

0*0 

0*0 

0*0 

0.0 

0.0 

0.0 

0.0 

4 


0.0 

0*0 

0*0 

0* 0 

0*0 

0*0 

0*0 

0.0 

0.0 

0.0 

5 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ROW 

/ COL 

1 1 

12 

13 

14 

15 






1 


73*4 

78* 1 

82*9 

07*7 

92.4 






2 


0.0 

0.0 

0.0 

0.0 

0.0 






3 


0.0 

0.0 

0.0 

0.0 

0.0 






4 


0*0 

0*0 

0.0 

0. 0 

0.0 






S 


0.0 

0.0 

0.0 

0.0 

0.0 










♦ T 

1 M E = 0.S000E402 S E C 0 

N D S 






AVERAGE HEAT 

transfer 

COEFF 

^ 




^ 





ROW 

/ COL 

1 










1 

0. 

S55E-02 










2 

0. 
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SHEET 3.3 INPUT/OUTPUT DATA FOR BLUNT-NOSED SOLID CONE PROBLEM 
(REFER TO FIG. 4) (SHEET 6 OF 6) 





Section 4 


SCRAMJET ENGINE SIDE-WALL GEOMETRY 


4.1 DISCUSSION 

This section describes the scramjet side-wall geometry which has been added 
to CAVE. An explanation of how the side-wall dimensions are used to construct the 
nodal network is given, followed by an outline of the input data format. 

A typical side-wall geometry is shown in Figure 7. The wall is assumed to be 
of unit depth and to be made of one material with constant properties . To facilitate 
generation of the nodal network, the side-wall is divided into three principal sections: 
the forward, center, and aft sections. Figure 7 shows how these divisions are defined, 
and their lengths are dimensions "A", "B", and "C" respectively. The user must 
input values for these lengths as well as the angles a and 0. 



1827-013B 

FIG. 7 TWO-DIMENSIONAL SCRAMJET SIDE-WALL ENGINE GEOMETRY 


The following discussion refers to Figure 8, which illustrates the arrangement 
of node elements in the scramjet side-wall as generated by Subroutine WALL2D. One 
should note that while the elements in the center section are all rectangular, some of 
the elements on the forward and aft sections are triangular in shape. Spatial incre- 
ments in the Y direction (AY’s) must be input by the user. The AY's are always con- 
stant in the horizontal X direction. Subroutine WALL2D uses the inputs for the 
dimensions "A" and "C," angles a and 0, and the AY's to compute the AX's for the 
forward and aft sections in the following manner. The triangular elements have one 
side forming a part of the wall's lower boundary. The angle of that side with respect 
to the horizontal is either a or 0. Because the AY for each element is a known input. 
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M ELEMENTS 



1827-014B 


FIG. 8 Nodal Network For Soramjet Side-Wall 



























































all triangular elements of the wall are right triangles , and one of the acute angles is 
known; the value for AX is determined. TheAX's in the center section are defined by 
the user, and are required input data. 

The number of spacial increments in the X direction (AX’s) is related to the 
number of AY*s the user has specified. The number of AX’s on the forward and aft 
sections of the wall are both equal to the number of AY’s chosen. Thus the total 
number of AX’s can be calculated as follows: 


Total Number of AX’s = (No. of AX’s in Center Section) + 2 x (No. of AY’s) 

Subroutine WALL2D performs a check on the wall dimensions to ensure that the 
sum of the individual elements truly represents the desired geometry. A discrepancy 
will cause an error message to be printed and execution of the program will be 
halted. Processing will be halted if the following relationships are not satisfied. 


Total Length of 
Scramjet Wall 
in the X 
Direction 


= A + B + C 
(Input) 

NDXMID 

i = 1 (Cente r Section) 

(Input) 

NDY NDY 

E AX 

i = 1 (Forward Section) i = 1 (Aft Section) 


and: 


(Computed by WALL2D) 

NDY 

A tan (a) = E ^ 

i = 1 ^ 

where NDXMID = Number of AX’s in the Center Section 

NDY = Number of AY’s used through the material 


The node numbering system is as shown in Figure 8. Nodes which appear to be 
located below the wall’s lower boundaries in the forward and aft sections are not used 
in the heat transfer calculations; they are used as an artifice for bookkeeping the 
node numbers. 
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For some areas of the scramjet wall, the user may want small nodal elements, 
such as the regions near the upper left and right corners. Because the AX’s in these 
regions are determined by AY and a or /S, the user must choose an appropriate AY to 
get a AX of the right magnitude. 

CAVE will permit the user to input a table of convective heat transfer coefficients 
and adiabatic wall temperatures for the top boundary of the scramjet wall. 

Alternatively, the user can let CAVE compute these parameters (Subroutine FLATH 
is used if it were computing H and for a flat plate) . If CAVE is to compute h and 

for the upper wall boundary, trajectory tables are required. In either case, h 
and for the upper boundary, h and for the lower boundary will be read in 
as Tables 2 and 3, respectively. When upper and lower boundaries are to be read 
in as tables, the upper boundary's h and T^^ will be in Tables 2 and 3, while the 
lower boundaries* h and will be in Tables 3 and 5, respectively. 

Problems with increased heating due to local interference heating or other 
effects during part of the trajectory can be handled using the HMODI flags. This 
option allows modification of the convective heat transfer coefficients along the 
upper boundary of the wall, and requires the user to input an additional table of 
heating modifiers. 

4.2 INPUT DATA FORMAT FOR SCRAMJET SIDE-WALL GEOMETRY 
Table Output Card 

• ITPRNT (II) 

ITPRNT = 0: Standard data tables will not be printed 
ITPRNT = 1: Standard data tables will be printed at beginning 
of printed output. 

Standard Tables Cards 

Standard Tables are read in so that the program can perform table look-ups 
of temperature vs enthalpy and pitot pressure. This method is used to obtain TTOT 
at Mach numbers >3.0 (thereby accounting for real gas effects). 

The input format for this table is described in the section dealing with Sub- 
routine NUREDl, given in Appendix A, Sheet A-1 of Reference 1. 
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• Table 1 contains: (46 cards) 

TTOT (°R) as a function (Arg. 1) and Log^^^ (P^ ^ (Arg. 2), 

\ o / 1 j 2 * 

where P. „ is the pitot pressure (in atmospheres) 
t, z 

• Blank Card (Signifies end of standard tables) . 

This table is provided as part of the computer deck and must be read in as part of the 
input data package whether it will be used or not. 


Indexes Card 




JGEO, L, M, NE 

JGEO 

L 

M 

NE 


3 (selects Scramjet Side-wall Geometry) 

Number of elements through the material 
(may be an even integer or an odd integer) 

Number of elements along the top surface of the 
scramjet side-wall 

Number of dominant eigenvalues to be used in solution 
(e.g. , a typical number is 5). 


(415) 


Geometry Cards 




A, 

A 

B 

C 

a 

0 


B, C, ALFA, BETA (refer to Figure 7) (8F10.0) 

= Length of forward section of wall 
= Length of center section of wall 
= Length of aft section of wall 

= Ar^le of wall boundary at front of wall (in degrees) 

= Angle of wall boundary at rear of wall (in degrees) 


• NDXMID 

NDXMID = Number of AX’s in the center section of wall 

• DX(1), DX(2), DX(3), ... DX(NDXMID) (8F10.0) 

DX(I) = Spatial increments in X direction for center 

section of wall (region of length = B) , ft 
(1=1, 2, ... NDXMID) 


• NDY 

NDY = Number of AY*s in the center section 


(15) 
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II IH nil 


• DELY(l), DELY(2), DELY(3), ... DELY(NDY) (8F10.0) 

DELY(I) = Spatial increments in Y direction for center 

section of wall, ft (I = 1, 2, ... NDY). 

Title Card 

• Run identification, comments, etc. (8A10) 

Radiation Card 


• EPSl, TBGl (2F10.5) 

EPSl = emissivity of surface 

TBGl = background radiation temperature, R. 

Material Properties Cards 

• MAT (15) 

• NMATl, RHOl, CONAVl, CPAVl (I10,3F10.5) 

• TMATl, TMAT1(2), .... TMATl (NMATl) ] omit (8E10.0) 

• CONDTl(l), CONDTl(2) CONDTl (NMATl) > if (8E10.0) 

• CPTl(l), CPT1(2) C TP 1 (NMATl) 1 nMAT1=0 (8E10.0) 

MAT = number of materials = 1 

NMATl = number of entires in properties table (maximum 

of 10). NMATl = 0 for constant properties 
RHOl = density of material 1, Ib/cu-ft 

CONAVl = average thermal conductivity of material 1 
(used when NMATl = 0), Btu/ft-sec-R 
CPAVl = average specific heat of material 1 (used when 

NMATl = 0), Btu/lbm-R 

TMATl(I) = temperatures in thermal properties table for 
which CONDTl (I) and CPTl (I) are given; 

1=1, 2 NMATl, R 

CONDTl (1) = thermal conductivity of materials 1 at temperature 

TMATl (I), Btu/ft-sec-R 

CPTl (I) = specific heat of material 1 at temperature 

TMATl (I), Btu/lbm-R 
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Initial Temperature Cards 


• KOBE, I, T(]), n, JJ (215, ElO. 0,215) 

• ... 


• • • 

• • • 

11100 


KOBE 

= 0 or blank 

I 

= node number 

m 

= node initial temperature, R 

the node number is incremented by the spacing JJ until 

II and JJ 

= the limited II is reached. Each node so specified is 
assigned the same temperature. 


Special Output File Card(s) 


• NPLOTS, ITCODE 

• NODNUM(l) , NODNUM(2) , . . . NODNUM (NPLOTS) 

(Omit if NPLOTS = 0) 

NPLOTS = number of nodes whose data will be written 

into special output file 

ITCODE = 1 if node temperatures in special file are to be in R 

= 2 if node temperatures in special file are to be in K 


NODNUM(I) 


= node number of temperature to be stored 


1=1, 2, ... NPLOTS. 


(214) 

(2014) 


Boundary Condition Cards 

Two options exist: 1) in the first option, the user inputs the flight trajectory and 
the code calculates the convective boundary conditions along the top surface of the 
panel in accordance with the equations presented in Appendix B of Reference 1, and 
2) in the second option, the user inputs directly the convective heat transfer coefficient 
and adiabatic temperature as functions of time and distance. 

OPTION 1. FLIGHT TRAJECTORY SPECIFIED 

• REFX, CODEX, HMODI (3E10.0) 

• GAM, RGAS, PR (3E10.0) 
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• NTRAJ (110) 

• TIMTAB(l) , TIMTAB(2) , TIMTAB (NTRAJ) (8E10.0) 

• ALTTAB(l), ALTTAB(2) ALTTAB (NTRAJ) (8E10.0) 

• VELTAB(l), VELTAB(2), VELTAB (NTRAJ) (8E10.0) 

• ALPTAB(l), ALPTAB(2) ALPTAB (NTRAJ) (8E10.0) 


REFX = effective boundary layer length, e.g, distance 

from leading edge of scramjet wall 

CODEX = 0. for uniform convective coefficient across top surface; 

CODEX = 1. for nonuniform convective coefficient 


HMODI 

GAM 

RGAS 

PR 

NTRAJ 
TIMTAB(I) 
ALTTAB(I) 
VELTAB (I) 
ALPTAB(I) 


(i. e. , function of x); CODEX = -1. for tabular input 
for convective coefficient and adiabatic wall temperature 
= nonzero value indicates that a table will be read at 
the end and used to multiply the convective coefficients 
= ratio of specific heats of air 
= gas constant for air, ft-lbf/lbm-°R 
= Prandtl number of air 

= number of points in trajectory table (2 < NTRAJ < 50) 

= time in trajectory table 1=1, NTRAJ, sec 
= altitude corresponding to time TIMTAB(I) , ft 
= velocity corresponding to time TIMTAB(I) , ft/sec 
= angle of attack corresponding to time TIMTAB(I) , 
in degrees. 


Convective Coefficient & Adiabatic Wall Temperature Tables 

Two tables are required to define the convective coefficient and adiabatic wall 
temperature along the lower boundaries of the scramjet wall. The first table gives the 
convective coefficient as a function of time (argument 1) and distance (argument 2) and 
the second gives the adiabatic wall temperature. The same values must be used for 
time and distance in both tables. The range of the distance argument must be from 
the upper left hand corner of the wall to the upper right hand corner as measured 
along the lower wall boundary, as shown in Figure 9. A description of the rable 
format can be found in the section on Subroutine NURED in Appendix D of Reference 1. 

• Blank Card (signifies end of tables). 
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I 




FIG. 9 LOWER BOUNDARY LENGTH FOR SCRAMJET WALL 


OPTION 2. CONVECTIVE COEFFICIENT & ADIABATIC WALL 
TEMP ERATURE SPECIFIED FOR UPPER & LOWER 
BOUNDARIES OF SCRAMJET SIDE -WALL 

• REFX, CONDEX, HMODI (3E10.0) 

REFX = effective boundary layer length, e.g. , distance 

from leading edge of side-wall, ft 

CODEX = -1. (indicates to code that Option 2 is being exercised) 

HMODI = nonzero value indicates that a table will be read at 

the end and used to mulitply the convective coefficients. 

Convective Coefficient & Adiabatic Wall Temperature Tables 

Four tables are required. The first gives the convective coefficient as a 
function of time (argument 1) and distance (argument 2), for the wall upper boundary. 
The second table gives the adiabatic wall temperature as a function of time (argument 
1) and distance (argument 2), for the wall upper boundary. In setting up these tables, 
the same values for time and distance must be used in both tables. The range of the 
distance argument must include the interval REFX to REFX plus the total length of 
the scramjet wall (i. e. , A + B + C). The specifics on the format for the tables are 
given in the descriptions of subroutine NUREDl in Appendix D of Reference 1. 

Similarly, the third and fourth tables should contain convective coefficients and 
adiabatic wall temperatures for the lower boundary of the wall. Distance must be 
measured along the lower boundary as shown in Figure 9. The range of the distance 

argument for the lower boundary must be from zero to ^ cQg (g) + B + 

• Blank Card (signifies end of tables) . 


COS (J3) 


)• 
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Time Intervals Cards 


• NTIMES (110) 

• NTIMES(l), NTIMES(2), TIMES (NTIMES) (8E10.0) 

NTIMES(]) = number of points in time intervals array 

(2 < NTIMES < 50) 

TIMES(I) = initial time (usually equals 0.) sec 
TIMES(I) = time at which temperature is calculated and printed 
1=2,3, ... NTIMES, sec. 

Convective Coefficient Modification Table (omit when HMODI =0.) 

One table is required to modify the convective coefficient on the aerodynamic ally 
heated surface. The multiplicative factor is given as a function of time (argument 1) 
and distance (argument 2). The writeup for subroutine NUREDl (Appendix D of 
Reference 1) gives the specifics on the format requirements. Follow this table with 
a blank card; omit the blank card if the table is not read in. 

4.3 SAMPLE PROBLEM FOR SCRAMJET SIDE-WALL GEOMETRY 

This subsection contains an illustration of the scramjet side-wall geometry. 
Figure 8 shows the nodal network. The section has an overall length of 11 ft and a 
maximum thickness of 0.5 ft. The boundary conditions on the upper wall are cal- 
culated by the code using a constant flight profile of 40,000 ft altitude and 3,000 ft/sec 

velocity. On the bottom surface the convective heat transfer coefficient was input 

2 

as 0.00555 Btu/sec/ft /R and the adiabatic wall temperature as 1,500 R. Both 
values were input to be constant for the entire time period of interest (200 sec) . 

The following pages show listings of the input and output data for this problem. 
The sequence is: 

Print -Out of Input Data: 

• Standard Table 

• Problem related input data. 

Print -Out of Output Data: 

• Geometry data 

• AX*s and AY*s 

• Volumes of nodes 
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• statement regarding storage allocated for S array in main program 

• Comment card 

• Material properties 

• Boundary conditions tables 

• Node number location within output array 

• Material number assigned to each node. (Only one material can be used with 
this geometry.) 

• The capacitance of each node 

• The conductance in the x-direction between nodes 

• The conductance in the y-direction between nodes 

• Initial temperature distribution. 

And then the following information is printed for each time interval: 

• Flight trajectory parameters, Mach number, altitude, velocity and angle of 
attack for each time interval 

• Average heat transfer coefficients calculated using the h values at the 
beginning and end of this time interval 

• Average heat transfer couplings, which include radiation effects, if any, 
calculated using the temperatures at the beginning of this time interval 

• Average adiabatic wall temperatures for this time interval 

• Temperatures at the end of the interval 

• Steady-state temperatures for the heat transfer couplings and adiabatic 
wall temperatures of this interval 

• Integrated heat input to each node. This gives the net heat transfer at each 
boundary node up to the end of this time interval. 

(Annotation has been added to the input and output to aid the reader.) 
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1890 


1 '0 ■ 

, ”51 ■ t ' . 

L 170. 


12/.0. 


1350. 1440- 153-3. 

l8-)0. 

1890 


1 • •.'•‘1 

■)■' » > j i.>i 'J t . 

U70. 


1260. 


1350, 1440. 1530, 

1000. 

109') 

•1 , • 

• 

!‘ji 18 

117./, 


1 260, 


135 ). 1 440, 153-3. 

1 80 Oo 

I89f) 

•1 . 

1'. ^ 1 

^ J'l'.. 

1170. 


1 260, 


135). 1440. 153->. 

10«)<». 

1890 


' ■; 


L L 7i). 


1 260. 


1351, 1440, 1530. 

1 799, 

301889 


. ~ 1 II 

-i ll 08 K 

L 1 70, 


1 260, 


135), 1440. 1529, 

1 799. 

1889 

1 ' 

' ■ 

■>.861 i8 ». 

L 1 7.). 


U6 0. 


1350. 1440. 1520-751790- 

761889, 

1 

1 • < 1 . 

"j 1 1 i 1 

U7o. 


1 259. 


135 ), 1433-581527-5 51 798, 

15108 7, 


1 1 

•‘1 f). -6 

2J.5 2.) 

.3 5,8(11 

2S. t 33.905 37,39 

40.2 

42.5 

- 

; 1 

: jt 

1 15-J. 

55 

324 J. 


36'-). 396‘J- 4105. 

4310. 

440-). 


1 . : * ' 

' ' • 2 6 2 1 . 

2881- 

12 

3242, 

093609.454003. 784210. 

4365, 

44 7.J. 

-1 

:>.• 

• •. >63 >, 


19 

3 2 0, 

613616. O7‘.O10, 844260. 

4440. 

4553, 

-1 . 

1 ■ 1 

• .. J J > . 

2 10 1. 

55 

3 2 '* 4 - 

233620-324 »36.65430<). 

4495, 

4600. 

- 1 . 

* 

. i.’. i . 

>8 1. 

8<» 

'24 5, 

2 

3622.684047. 144320, 

4505. 

334637, 

- ‘ 

1 

J 'lJ ). 

MB 1. 

8 

32‘i5„ 

7 5 3625. 5 14 J62.6 743 50. 754 551. 

234695, 


,ti .}r',b.Z5i6Z r,ij7'»')r6.il843 79.^ 7-« 5 95. 6 7 54 . 55 23 T TO T 

.3»i3:4o, 2 73fa?'7.294.>8 5, 194394. 524 625. 084 79 3. 7 224 T TOT 
.90 3 246. 2 73629, 764088.5 44 40 7, 7 1464 0. 5448 1 5. 1 92ST TO T 
.8 I 3246. 2 7363 J. 714 >94. 9 244 22 . 894 665. 7448 45. 9526T TOT 
o35j 246.2 7363 I, 184099, 9M44 34, 92468 7. 384880. 49 2 7T TOT 
, i:i2*6.2536M. I 6 4 1 01 . 62444 t . 754 70 1 . 2 7490 1 . 2 Z 267 TOT 
.^8 3 245. 7 536 3 >.864 102.8 54445.994 70 7, 7449 10. 7 729T TO T 
,09 3 244,2 5362 9.864 10 3. 2 844 4 9.914 71 7.074926.1 1 30T TOT 


STANDARD TABLE 


5 1,) 

() ), ) 

7 1,0 

B5, ) 

9.).8 

^ -j / 1 . 

/ J-i5. 

5 1 7.‘- 

5.35. 

5 75 0 


• -65, 

4JJ0, 

55 75. 

5925 

» -t 7 • , 

■*) 7 5, 

5 i(.5. 

5745. 

6 )65 

■» # t ^ . 

>29 .*• 

557 >. 

5 865- 

6185 

1/15, 

3 34 

5645, 

5)30. 

6 260 

j »■> 'i 

>i5 5. 

5;>JJ. 

6 11). 

6465 

5 1 i } m 

.> J J5. 

554 5, 

6265, 

6615 

5 M * . 

.7 )•'. 

o >6 5, 

641 ). 

6 74 5, 

52‘!'i.. 

i 740, 

6U5o 

6515. 

6«8'J, 

J ) J ) • 

J 18 5, 

6 310. 

66 73. 

7-I5-I. 

^ 1 2 I » 

>795. 

6468, 

685 ). 

7255, 

5 '♦ 3 , 

:• 1)5. 

658-J, 

7313. 

7425, 

J i * 5 • 

> 1 1 3, 

u ^ 'i « .’ • 

7 >85. 

7515, 

^ 32 ' » 

.2 '■), 


7 26 ). 

9599' 


J^E^LyMyNE 


,0 0 ;. ) 9.>,0 10 3,0 IIO.O 120,0 31TT0T 

7.'. 5 .35. 5 7S'J. 6 3 25. 6905. 755 J. 321107 

JO, 55 75. 5925, 64J0. 7120. 7665. 33TT0T 

55, 5745 . 6)65. 6530, 7195. 99999+934T TOT 

73 . 5865- 6185- 6610, 7260. 99999*-93 5T TOT 

45 . 5330, 6260. 6690, 7305. 99999*936TTU T 

JJ. 61 1 ). 6465. 0845. 7435. 99999»-93 7T TOT 

45 , 6265, 6615. 7025, 756'l, 99999f 930T TOT 

;o- 6411. 6 745, 7180, 76871. 99999*-939TT0T 

45o 6515. 6880. 7295. 99999f999999*940TT0T 

10. 66 73. 7-l5'l, 7480. 99999»-999999*94lTT0T 

6il« 685 ). 7255. 7580. 99999*999999*942TT0 T 

a-J, 731 3. 7425. 99999 + 999999*999999*943TT0T 

3.'. 7 185. 751 5. yV999f 999999 ♦999999«-944T TO T 

7 26 ). 999 99»9<;9999f99g999 +9999 39*945TT0T 

BLANK CARD 

-. 46 2 3 14. 036-^— -A. B. C. ALFA. BETA 

NUMBER OF AX'S IN CENTER SECTION 

'.5 .1.5 ».5 0.5 AV'g 

NUMBER OF AY'S 

. 1 O.l 0.05 ^Y'S 


1827-016B 


SHEET 4.3 INPUT/OUTPUT DATA FOR SCRAMJET ENGINE SIDE WALL PROBLEM (SHEET 1 OF 11) 





SCRAHJFT wall v^EO.Mt^RY TEST C iSfc WITH Iqi* oOUNOiRY m 

0> ■), ^ - .,... 

I ^ 

lo^* \j, 2 IS< ^ — 


I 'iJ.i.'J 

I 1 1«)0 
ll I 

1 3 ll 2<t 2y 

' Oa 2 s J '} a 

1.4 33.3 

i* 

• 1. 3. 

<«.tJJ3a ^JJOlJa 

3)10. 3rJi). 


4 r d2 <JJ 95 96 




M. JJ« 

W'J;).) . 
3iii)0. 30-I I. 


i.Mii MW cvicDi COMMENT CARD 

MAT 

INITIAL TEMPERATURE 

SPECIAL OUTPUT FILE CARDS 

AIR properties 

•FLIGHT TRAJECTORY 


EPSI.TBGI 

MATERIAL PROPERTIES 
SIGNIFIES END OF INIT TEMP CARDS 
REFX, CODEX, HMDDI 


« 0» J 

12-32 sew 4 ^ JET WALL INAfh' SLRFACP hT ')IST. 

2JJ. 

#•1 a«Jil?5‘j .>]iJ3‘33 

.3 .00535 . 03355 

. 3.2 .00555 .JJ555 

•203 TAW SCKAMJRT IM'^ER WALL SURFACE OIST. 

U..) 2»j« 

.■J 1500- 15)). 

.,•)') 15')0. 15 0 >. 

-).2 l5jJo 1 5 JO. 


LOWER WALL BOUNDARY 

HAND Taw distribution 


-TIME STEP INTERVALS 


-BLANK CARD 


1827-017B 






JLfA- 

*Lrrt= .jtC-itES 

t_M= 

KclAs 'C,(>iL£- 


PRINT -OUT OF INPUT DATA DESCRIBING SCRAMJETWALL GEOMETRY 


-LIST OF AX'S IN CENTER SECTION OF WALL 


0.50 0«=0 


hOH'dCJ- :■= = 


-LIST OF AY'S 


CtL«»S FC-; , 

0« 3( J 0«o0 J 

0 • bO J 0 • 500 

0*400 0*<iU0 


CELY'5 FCfi *»t-L 


LIST OF ALL AX'S - INCLUDES THOSE INPUT FOR 

CENTER SECTION & THOSE COMPUTED BY CAVE 

O.jOJ C.500 0*500 0*500 0*500 

0.500 0.200 0*400 0*400 0*400 

LIST OF AY'S USED IN PROGRAM (SHOULD BE SAME 

AS ABOVE) 


1827-018B 


SHEET 4.3 INPUT/OUTPUT DATA FOR SCRAMJET ENGINE SIDE-WALL PROBLEM (SHEET 3 OF 11) 




K/O.L 


J*JJJ JaJcC 


. w . : ^ 

• J 2 J J • J 


OaJ>^3 JaJcs JaJJt; 0*025 0.0 25 0.025 0.025 0.025 0.025 0. 025 

J.OSO 0.05o 0.050 0.050 0.050 C.C50 0.050 0.050 0. 050 

0, 050 0.05 0 0.0 =o C.05J 0.050 O.OSJ J.JSO 0.050 0.030 0. 050 

:.w50 J.u5. J.us: 0»0 5 0 0.o:.0 0.050 0.050 0.05C O.OSO 0.050 

^.uIjO O.wiO 0.0 = 0 J.050 0.050 0.050 0.050 0.050 0.050 O.OSO 

O.O^t 0.02s 0.0^5 0.025 0.025 0.02b 0.025 0.025 


VOLUMES OF NODAL ELEMENTS 
COMPUTED BY CAVE 

NO OF ELEMENTS ON SCRAMJET 
/WALL SURFACE BOUNDARIES 


NO. OF NODES 


vNO.OF EIGENVALUES 

3 (^OUH INANT MOOES...REQUIPES 6644 WORD 

%F.05 ARRAY SIZE REQUIRED 

FOR THIS PROBLEM 

t **^*******m*■*tmt**■**■¥*%*^*^‘*^^**m*■^■^**■»■■^■^■*■ *■*%**■»■* 


6644 WORDS OF MEMORY 


SHEET 4.3 INPUT/OUTPUT DATA FOR SCRAMJET ENGINE SIDE-WALL PROBLEM (SHEET 4 OF 11) 




SC«iMJt7 4 *lL CEC^clt', itST CA:t -Mr- luP CCvtCA^-V 


lAk call ^L. A TEOa 


•» • k. j Lti .‘ / C v-»- T 


LP=0.i0uw ylU/L G .« 


CONSTANT MATERIAL 
PROPERTIES 


Tl»<t IN CcCCNJG 

Lfm 'jJJsiC *■ 

alt ITO.^t iN FcdT 

S. 4 J 0 JL»C~ C-.-ujCttC' 

VElCCITV IK F=u( Pr.i^ otC» 

J »JCJJP •C-4 J»3J 
AKUL, CP JiTACiv In 


TRAJECTORY TABLES FOR UPPER 
SCRAMJETWALL BOUNDARY 


SCHAHJLT *All iNKCi! SLSfA^; 


r*m FCK cCTIC*' SUfcF ACC 
Hr DIST. 


H& Taw tables READ 
\H FOR SCRAMJETWALL 
LOWER BOUNDARIES 


0.0 

0*0 

..2CCCEACJ 

0.0 

C • 'J 

0.0 

0.0 

0 .0 

OaO 

0.0 

v).0 

Oa bbLiOc- 0£ 

CaStEOE-Ci: 

Oa 0 

c.o 

OaO 

OaO 

OaO 

OaO 

0. 0 

0««Oa3E«0 1 

0«SbtJc'J2 

Ca5E50E~02 

Oa 0 

0.0 

OaO 

OaO 

OaO 

OaO 

0.0 

0* tosoE+Oi; 
2 3 

0 .EEEOE-02 

0 .&SS0E-02 
TA» SCRAt'j 

O.C 

ET INKER 

0 .0 

kALL surface 

0.0 

OIST. 

0.0 

OaO 

0.0 

Oa 0 

0.0 

il.O 

Ca20C0E*03 

O.C 

0.0 

OaO 

OaO 

o;o 

OaO 

0.0 

J>0 

0 . l!iiJuE«-0« 

C. I5C0EAC4 

0.0 

0.0 

0.0 

0.0 

OaO 

0.0 

OaO 

«>•* COOE 1 

!)• I!>00e«04 

Ca I5C0EAC4 

OaO 

OaO 

0«D 

OaO 

OaO 

OaO 

OaO 

0.1 020EA02 
0 0 

0* ISOOEFOA 

Oa ISO0EAO4 

0.9 

o.o 

OaO i 

OaO 

0.0 

OaO 

0.0 


1827 - 020 B 


SHEET 4.3 INPUT/OUTPUT DATA FOR SCRAMJET ENGINE SIDE-WALL PROBLEM (SHEET 5 OF 11) 



NUVcEPS 


NODE ARRANGEMENT - RESEMBLES 
SCRAMJET WALL GEOMETRY 


RCK / CCL 

1 

_ a 

J 

.. 

5 

c 

1 


t,V 

1 i. : 

1 o 

2 5a i) 

j 1 aO 

2 



LaO 

14.0 

20. 0 

£ O a 0 

32.0 


OaO 

2 a 3 “ ^ 

1 E. 0 

2 1.0 

27. J 

Jj.O 

•* 

0. J 

C ,0 

C • u ' ' 

.-‘C. y 

2da 0 

J4a0 

5 

OaO 

C. J 

C.O 

o7^ 

2>. J 

35.0 

6 

0. J 

0.0 

0. J 

J. j 

oTir 

Jfc. 0 

RCa / ecu 

1 1 . 

12 1 

3 

1 4 

13 

1 c 

1 

C la J * 

o7.0 

73a J 

7S.0 

bS.O 

91.0 

2 

c2.0 

EE.3 

74.0 

*10. c 

9t a 0 

92.0 

3 

l3. J 

CS .C 

7 E* 0 

0 

L 7 . y 

93.0 

(• 

C4a J 

7 J. J 

7c. J 

ac • j 

' ■ « J 

94a 0 

E 

O5a0 

M.u 

7 7.0 

U3 J 

f'» a J 

95.0 

6 

CO a J 

7a .0 

•' £ a_C 

S'* J 

_ 90a_J 

96.0 

RCk / CCL 

21 

^ 2 





1 

12 1.0 







2 

122.0 - — 

o.c 





3 

— - — X)«i3 

OaO 





4 

0.0 

0.0 





5 

0.0 

c.o 





« 

OaO 

OaO 





waTLAImu 


NCtiL 





PC'* / CCL 

1 

a 

J 


- 

3 

1 

1 a J 

1 aU 

1 . 0 

L a J 

1 .L 

1 ... 

2 

J « V 

. 

1 . 3 

. • - 


! a - 

2 

0 a J 

LaO 

1 • o 

1 a 0 

1. J 

U J 

4 

\Jm J 

U. J 

0. 0 

1 a 0 

1.0 

1 .0 

5 

Oa ^ 

C.O 

J. J 

0.0 

1.0 

1 aO 

C; 

0.0 

c.o 

C. J 

OaO 

.'a J 

1.0 

RC'A / 

1 1 

12 

13 

14 

15 

ic 

1 

l.J 

1 .0 

1 . J 

1 . J 

1 .0 

1 .0 


. > 1 

1.0 

1. C 

1 a J 

t .0 

1 . J 


1 . J 

1 .0 

1 . J 

1 • 0 

1.0 

1.0 

4 

1 . J 

1 

1 . c 

1 • y 

1 .0 

1 .0 

3 

1 aO 

1 au 

1 . c 

1 a 0 

1.0 

la J 

tJ 

i av' 

1 .0 

l.O 

laO 

1 aO 

l.C 

fiCa / 

1 






1 

laO 

. aO 





2 

1. J 

C.O 





3 

0.3 

C.O 




1 

4 

J lO 

CaO 





5 

OaJ 

C.O 





t 

OaO 

CaO 






_7 

e 

^ 

1 3 

37. 0 

43.0 

4 9.0 

55. oT 

3Ba 0 

44a0 

50.0 

56.0 

39.0 

45.0 

51.0 

-■7. J 

4U.0 

46.0 

52.0 

56.0 

41.0 

47aO 

53.0 

59.0 

42.0 

__4e.o 

54.0 

60.3 

r 7 

1£ 

_19 

20 - 

97.0 

103.0 

109. -3 

llSaO 

9d.O 

104.0 

1 10. 0 

116.0 

99.0 

lOS.O 

1 1 la 0 

1 17.tl- 

tOOaO 

106.0 



— 0.0 

101.0 

107.0 

0.0 

0.3 

102.0 

— fflo 

0. 0 

0.0 


ua 

la 0 
lai) 
la J 


1 a J 
laO 
la 0 
laO 
laO 


1 a i 

1 #a 

laO 

I a 0 


la J 
laO 
la 0 
laO 
laO 
OaO 


laO 

laJ 

i.a 

laO 

laO 

laO 


laO 

laO 
1.0 
1 . 0 
0.0 
0.0 


ONE MATERIAL THROUGHOUT 


1.0 

laO 

1.0 

1.0 


laO 
la 0 

uo 

OaO 

OaO 

OaO 


1827-021B 





Oi 


OhACITASCc LACt~ Ni;je 


RCM / LCU 

1 

< 

J 

« 

5 

t 

7 

e 

9 

10 

t 


a« IJnJ 1 

Ja (CLr: 'bt 

b J 1 1 i.- -1» . 

... i JbL • j; 

Im is. ,'E *\JJ 

J.« 74f;»J J 

J. 0 7<;r^ua 

jad7<;E*aa 

Jaa7«c*0d 

J 

J *0 

w, l«br-* -l 

J a c 1 ' II 

-a a. 1 lE *V- » 

J. bl IL * - 1 

J a 1 J ^ t' 1 J 1 

4a t 7 uE au 1 

J. 1 /6t *01 

Oa 1 76E*0 l 

Ja 1 76e«-tf 1 


0 • b 

w . C 

J . i « ■ ► J 1 

b * c 1 1 1^ « ^ 1 

J a b 1 1 h • ^ 1 

J a ! J jE *»/ 1 

3a 1 7tc f V 1 

Oa 1 7fcE*OI 

Oa t 7<>£ *ii 1 

Ua 4 7«c*a I 

4 


J • 0 

imsl 

b# 1 t bK *-u 1 

J. 1 iL • J 1 

Ja . 

Ja i rbE«^Jl 

0. 1 76E *0i 

3al7be*0l 

Ja 4 7b£* 01 

b 

Q.u 

3. J 

..b 


■ Jm . JbL * ■.' 1 

J a 1 0 1: • U 1 

b. 1 7&E*ai 

1 7bE «0 1 

3a 1 7f;b »0 1 

J. 4 76E*dl 

t- 

aO 

4 • w 

: » .. 

b . ■: 

.. . - 

a • a. O b *0 <j 

b ad7'?E4>00 

0. d 7 <;e *^00 

Jad79E«^00 

J. d79E* 00 

RO* / CwL 

1 1 

1 c 

J 

1 4 

1 b 

1 c 

1 7 

le 

1 9 

23 

1 

O.C 7SL*V J 

jac7>E*0J 

b a E ' J 

Cat 7CE ♦ J J 


c ac 7<It »0 3 

Oa 30 

0. 7 J3E »00 

6a 70 3E *00 

0. 703E«-00 

2 

3 • 1 7t»h » 3 1 

Z , iTtE^O 1 

C. 1 

>* . 1 ? c E ♦ b 1 

Ja 1 7t. E » «) 1 

J . 1 7uF < J 1 

0 a 703t*33 

Oa M 1E*0 1 

Oa l« IE *0 1 

Oa 1 4IE*01 

3 

Jal7cb*01 

0. 1 7LE * C 1 

Ca 1 Itc 

ba 1 IcC 1 

Oa 1 /oi: • 0 1 

J. 1 7oF *0 1 

Ja 7 33F «v) J 

0. 1 4 IE *01 

<3. t 41E*0 1 

0a7C3E*03 

4 

> • 1 * J I 

b . 1 7lL * 0 1 

J . 1 :tL 

Ja 1 7Cl- iJ 1 

Ua 1 7ut 1 

.■ a 1 7c e • C 1 

J. 7 J3E03 

0. 14 1E*0 1 

Oa 70JE*00 

0. J 

5 

C . 1 7at * 3 1 

J* l7Cii*0l 

J . l 7cE ♦ J 1 

V: . 1 7 »■ c ♦ 0 1 

J, 1 7nC *0 1 

Ja 1 /oE^b^l 

J. 7 JJE*33 

Oa 7J3E400 

Ja 0 

OaO 

t 

Jfc *U J 

J, ^7St J 



-■a ' 7 ‘at a J J 

1, t. *.i J 

Oa 1 7uc » J J 

Oa'O 

Oa J 

<3.0 

HC« / CL.L 

. 1 

w • J -1. • J J 

• 1 / c e * J 


\ZER0 VALUES, 

BECAUSE THERE ARE NO NODES 




m 7 bob 
C • J 

wO 



AT THESE LOCATIONS 






C *U 

1 1 ^ ^ 

— ^ 








•j 


. ■ b 









C 

3, si 

j •«: 










CC ^CULT .SK El (S M-Ul 

.^t<. I IC N 









RCV 

I' 

t-OL 

1 

2 

J 

•4 

b 

6 

7 

e 

9 

10 

1 



0 a b - J J 

J . J t sE - 0 2 

J ■ JtsE-02 

b. JL5E- J2 

ba<*7tJt-J2 

0. = JJE -02 

0 .43 4E-02 

0.4JIE-02 

0a*3 lE-02 

0.431E-02 

2 



OaO 

3. 74 cE - 

Ja 71 cE - 32 

Ca 7 ldE-U2 

Ja <«b7E -u2 

J. 1 JJE-0 1 

Ua db lb-32 

Oa bb lE-02 

0.8blE-02 

3.dolE-02 

3 



b a L 

0 a *3 

Ca 7 1 tc-02 

wa 7 1 2E -02 

Oa sbre-oa 

Oa lOlt -0 4 

O.at lE-02 

O.eolE-02 

0a8btE-02 

OaablE-02 




Cab 

OaJ 

0. j 

0. 74tc-J2 

0. 5»S/r-0 2 

0 . IJ-JE- J 4 

0.86 4E-J2 

OaSblE-02 

OaBb lE-02 

Oa86IE-02 

S 



d a «i 

0 a C 

Ua 0 

Ca L 

Ca ‘.37 E- j2 

3. 4JbE-0 1 

O.BblE-J2 

0. 0b lE-02 

OaBblE-02 

0.e6lE-02 

c 



0.0 

b a br 

3. <3 

Oa 0 

Ca 0 

J. £j-3t- J2 

d.4j lE-02 

Oa«JlE-02 

Oa*31E-0Z 

0.4J1E-02 

nc« 

/ 

LCL 

1 4 


1 J 

1 4 

1 b 

Ic 

1 7 

18 

19 

20 

1 



Oa«l j It- Ja 

0.4 J 4.1-C2 

0.4 J l^-C2 

Ca 4C 1 1- 02 

J.-JlL-v . 

J a c. 1 jE -u 2 

Ja 710E-O2 

0ab3EE-J2 

0a53e£-02 

0aS3eE-02 

'2 



ead6lb-02 

3ac6 lc-02 

3 • 1 c - 02 

Oa eel t ^02 

Oa do] t-02 

Oa 12JE -i3 1 

Oa 144E-01 

Da 1 3aE-Dl 

Da 108E-0I 

Oa lOBE-Ol 

3 



OatJt lE-02 

0 a€c 1C-C2 

0 aCc it-ca 

C. Cc lE-Oal 

C.BolE-02 

0. 12JE-31 

Oa 144E-01 

Oa lOBE-OI 

Oa lOSE-DI 

^OaO 

4 



0aC61E-32 

OaEblE-02 

Oatblc-02 

Oa EME-02 

0. EblE-02 

0. 12JE-0 1 

0. 144E-01 

OalOEE-OI 

^.0 

^ OaO 

‘5 



Oaft6lE-02 

OaEblE-02 

0«ebtE>02 

O.et lE-02 

O.eo lE-02 

Oa 12JE-0 1 

Oa 144E-01 

aOaO y 

^0.0 

OaO 

6 



0a43IE-02 

Ca«3 lc-0« 

Ca4J tc-OZ 

Oa 4 2 1E-U2 

0.43lE-02 

, O.e 1 9E-02 

0 a J 

4o.o y 


OaO 

A0» 

/ 

COL 

21 

22 





w 

/ 



1 



Oa7ldE-32 

r 





A i 

f X 



2 



0.0 

c.o ' 









3 



C*0 

0.0 

- 

~ 



\ / 




4 



0.0 

OaO 



— 


\ / 




S 



0.0 

OaO 




— 





6 



OaO 

OaO 




' 

\L/ 












ZERO VALUES, BECAUSE THERE ARE NO 



NODES TO THE RtGHT OF THESE 
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^.L^L;uc^^^Cc I.s T-JIfi 
ACa ' Ct- 1 

c 1. T iL N 




C. 

7 

e 

9 

to 

i 

O.u 

J.cStttJJ 

J. 2£ce J 

i • 2 t t L * 2l Li 

Ja 25iuh *■ 2c; 

Oa Ijyfc+OC 

0.21 £E*00 

0.21£E*0C 

0.2:i 5E*00 

0.215E*00 

<■ 

0.0 

J • «- 

0.c::c:; *i;j 

1 

2a c'jfiE *00 

0. 1 2-7C ♦ Ju 

J.2l5E*00 

0.2I5E*00 

0.21SE*00 

0.215£*03 


0 • J 


.1 * ^ 

2. 2lc.i.*JJ 

0. 25dE 400 

la 1 2 a t. 4 0 lI 

J.216E*0 0 

0.21 £E*00 

0.21 5E*00 

0.215E+00 

4 

J • ^ 

J.O 

J. J 

2. 2 

2*i.Scc4i.J 

■>*» l^'.L*Oj 

0.2 1 6“ *03 

0.21 £E*00 

3.2i£E*03 

0.215E*00 

- 

J.O 

. : 

:■ . 0 


2.2 

C..1 /n::,40iy> 

Ja20 7S^*jO 

J. 28 75*00 

0.23 7E*00 

0.26 7£ *00 

c 

J d J 

- > - 



- . J 

W* • J 

0.0 

0.0 

0. 1 

0.0 

fii:« / tcu 

1 i 

12 

S J 


1 S 

1 ■- 

1 7 

IE 

19 

2 3 

i 

0 l,:i: Tv/vi 

- . ^ ; d L - . 0 


i jL * J ^ 

.. 2 1 3E 4 . w 

Jadl'J E*0o 

Oa 3ulE~0 I 

Oa 1 72E * 00 

0. 1 72E*0a 

0. 1 72E*00 

2 

li.b 1 3C. *-wv) 

J 

J. 2 1 * i. J 

2*21 Jii *■ LI 2 

J. 2 1 jf *Ow 

Ja21-JC*Uj 

0. 3olE- J 1 

0. 1 72E*00 

0. 1 72E*00 

0. 1 72C400 


J • 4. 1 JL! 4 y J 

J 


kiB4 l=l.r^t. 

C i r ^ U 

J»£ 

u a o ^ » t- - 1 

wal<*2E*uO 

0. 17 25 *00 

J.2 

•* 

J . 2 1 '-Jc ▼ 0 u 

J ■ 2 1 L E 4 C 0 

0. 4 lEE ♦ j: 

w«2 13E*02 

2. 21 se. »2 J 

Ob 2 laE.uO 

0, ac iE-0 1 

0. 1 72E*00 

0.0 

0.3 

£ 

tC /t 

fc*00 

J.4^ 

.*4.t ?E*02 

J. 2U7t *0.i 

J. Ed 7t*0 J 

0. 1 ISEfOJ 

0.0 

0.0 

0.3 

b 

J ...• 

J » ^ 

0.0 

J 

C« 0 

0.0 

0.0 

3.0 

0.0 

3.0 

RL« / CLL 
1 
2 

t 

IN^T I AL 
uk ^ ccl 

i L 

C« 1 *00 
C.O 
C .0 

c * c 

1 1 u A T o-'. t: 
1 

0 ■ V. 

J.L 

c ■ c 

lil ^ li;w1 ICN 
2 

Jt.C .4 


ZERO VALUES BECAUSE THERE ARE NO NODES BELOW THESE 

t o 7 e L w 

1 

^ JO. J 

£00. J 

5 0 J. 0 

bJ J. 0 

SOJaO 

500.0 

500.0 

500.0 

SOO.O 

500.0 


• J 

£u0« J 

£C0a C 

SCO. 0 

500.0 

500.0 

500.0 

500.0 

500.0 

500.0 

- 

0*0 

C*0 

£C0a C 

£ C 0. 0 

SOOa J 

500.0 

500a0 

£00.0 

SOO.O 

500.0 

4 

li • O 

<1.0 

0. J 

£0 0. 2 

5C0.O 

500.0 

500.0 

£00.0 

500.0 

530. 0 

£ 

W • W 

w • J 

c. : 

0. C 

sue *0 

500.0 

500 . 0 

£00.3 

SOO.O 

530. 0 

6 

J 

LaO 

Oa 0 

J. D 

0.0 

500. U 

SOO. 0 

500.0 

500. J 

500. 0 

Cm / CCt_ 

i 1 

12 

1 0 

14 

15 

1 (J 

1 7 

le 

19 

20 

1 


£C0. J 

52 J- 0 

5t 2. W 

SCO .o 

5 0 2.0 

500 ■ 0 

£00. 0 

500.0 

500. 0 

2 

jo 0*0 

£ C 0 a 0 

e . 0. 0 

SO 0. 0 

SOUaO 

500.0 

500.0 

500a0 

SOO.O 

500. 0 

2 

0 

di. Ja J 

: o . J 

£0 0. 0 

500.0 

500.0 

500.0 

500.0 

500.0 

500.0 

% 

buua 0 

£C Ja 0 

5C2. 2 

S L 0 a 0 

5C0.0 

5 0 0.0 

500. J 

503.0 

SOO.O 

0.0 

S 

UU 0. 0 

L- 0 C a J 

5C 0. 2 

bOOa U 

500.0 

500.0 

5 OOa 0 

£00.0 

0.0 

0.0 

6 

i-OOaO 

£00 a 

£00. 0 

SOO. 0 

500.0 

5U0.0 

500.0 

0.0 

0.0 

0.0 

a\i / ccL 

1 

2 

3 

4 

5 

6 

«>« * 1 

21 

oOOaO 
t> 0 J • 0 

0* J 

i).0 

0.0 

1 *•* + *• *****1 

£COaO 
C a 0 

CaO 

0*0 

Oad 

OaO 

k*«4«********h*****»4***-»*4*** + ***»***4*l-+4**4-** 



I- 4 • * f + « t * « • 4 
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ALMTCCb^ aCCOj.J 


VCLCCMts JCdiuOO 


angle CF at T ACK9 0»0 



tAl JMAN^Ft" 

CCtff 1C lEM 

s tju/sec- 








R0« 

^ CCL 

1 

i 

3 

ft 

5 

b 

7 

0 

9 

10 

1 


it 1 I7c*0l 

0«ltoE*0t 

0* M4E-0I 

0* 1 JOE-01 

0. I07E-0I 

0. lOJE-OI 

0. lOOE-OI 

0.57IE-O2 

0*0ft2£-02 

0* .fls£-u2 

2 


J.J 

o.*:ffsE-o2 

0.0 

0. 0 

0.0 

0.0 

0.0 

0.0 " 

0.0 

0*0 

3 


b.O 

o.c 

C.SSSE-C2 

c.c 

0.0 

0.0 

0.0 

0*0 

0.0 

0.3 

4 


c.o 

o.c 

0.0 

O.StbE-oa 

0.0 

0.0 

0.0 

0.0 

0*0 

0*0 

b 


0. J 

0.0 

0.0 

0.0 

0.5t»5t-02 

0.0 

0.0 

0.0 

0*0 

0.0 

b 


o.k. 

o.c 

0*0 

c.o 

0.0 

0.S!>bfc-02 

O.SSbE-02 

0.5SS6-02 

O.SS5E-02 

0.555C-02 

P0» 

/ CCL 

1 1 

12 

13 

1ft 

1S> 

1 1 

17 

IC 

19 

20 

1 


C.b>rf0t-02 

0.«ttd-02 

o'.04 4e-o2 

0. 92ftE-0? 

0.e03E-02 

0.7d7t-02 

0.771E-02 

0.766E-02 

0*742E-02 

0.729C-02 

2 


C.O 

o.c 

0.0 

0.0 

0.0 

O.C 

0.0 

0.0 

0.0 

0*0 

3 


0.0 

0. 1 

0.0 

c. c 

0.0 

0.0 

0*0 

0*0 

0*0 

0.555E-02 

4 


0.0 

O.C 

0.0 

0.0 

0.0 

0.0 

0.0 

0*0 

0*559E-02 

0.0 

s 


0.0 

o.c 

o.c 

0.0 

0.0 

0.0 

0.0 

0.S9SE-02 

0.0 

0*0 

« 

-W0» 

1 

2 

J 

4 
9 

5 

/ COL 

0*99SL~02 

0«717E^02 

0.955E-02 

0.0 

0*0 

0.0 

(T.O 

0.S53t-02 

££ 

0.705E-02 

0«0 

0.0 

0.0 

0.0 

O.C 

0.&5SE-02 

0.5S5C-02 

0.55Sfc-02 

0.55SE-02 

0.355E-02 

0*0 

0.0 

0*0 


AVtI-AOE “f 

-nNiPiin 

^ 1 NO 3 

LTL.'SEC-UbO 

. H 






rtUH / 

LCu 

1 

i 

J 

4 

b 

b 

7 

Q 

9 

1 


w* ^ ! Cc-u J 

C.cCCfc-02 

0.CU2E-02 

.. • ob3t-0<: 

J.oftlb — ‘.'2 

3.21 >c-02 

0.03IE-J2 

3.ftdCE-32 

3.471E-02 

2 


0*vt 

J*-JFE-02 

J* 0 

J. u 

3.3 


0.3 

3.3 

3. J 

3 


c..> 

0 • c 

0. J3feE-02 

c.c 

0.0 

0.3 

0*3 

0.0 

0*0 

ft 


0*3 

0*0 

J.O 

o.33ee-02 

0.0 

3.3 

0.3 

0.0 

0.3 

S 


0*0 

0 * c 

0.0 

J.O 

0.3J0fc-0C 

0.3 

0.0 

0.0 

0*0 

6 


0*0 

3*0 

0.3 

0. 0 

0.0 

0. lb •t-02 

0.277E-02 

0.277E-02 

0.277C-02 

ROft / 

CCL 

1 1 

IC 

13 

ift 

l!i 

to 

1 7 

le 

19 

t 


U.ftftbd- Jw 

J.-.1JL-C2 

3.422C-C2 

C. ft lftE-02 

0.4 J2I--32 

0.39,€-Ja 

0. lSftE-32 

3.302E-02 

0.297E-02 

a 


C*J 

. w 

0.0 

0.3 

0*0 

0*3 

0.3 

0.0 

0*0 

3 


C * w 

3 . \ 

J. C 

w • C 

0.0 


3.3 

0.0 

0.3 

ft 


3* J 

J« . 

3*3 

C. 2 

3. J 

0.3 

0*3 

0*3 

0.2C9C-02 

5 


0.0 

. • J 

3 . V 

J. 3 

3.3 

3.3 

0.3 

0.229E-02 

0.3 

c 


w . ft / 7u- J- 

; 7E- 

C.2?7£-C2 

C. 23 7r- J.. 

3* 2?7h-3C 

3...7!E-3J 

C. 1 lftE-02 

0.0 

0*3 

RC* / 

CCL 

<1 

ft ft 








1 


C . ft 2 - 

ja . 1 1. - 2 








2 



* 








3 


C.O 

3 . w 








ft 


C*0 

• 3 








« 


0*0 

C *0 








b 


C .0 

« . . 









to 

J.4SrE~02 

OmO 

0*0 

0*0 

0.3 

0«277E-02 

20 

0-29JE-02 

0«0 

0.2296-02 
0. 0 
0.0 
0.0 
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5TF Ale 

TLmPEhaI v/he.7 

FCh IFF 

BljLNCApy CLNCITICNS Af 

THIS time 

uto.n 





C« t Clu 

\ 

2 

3 

•4 

b 

b 

7 

B 

9 

kO 

1 

lie 7,ti 

1 1 7' .2 

1 1 { ;• n 

1 1 ^0* b 

1 1 93.9 

1 1 9e* 1 

1 19a. 2 

1200.8 

1203.6 

1206.4 


J • J 

1 1 /<;• o 

i 1 ds.** 

I 1 (>4. « 

1197.7 

1199.9 

1201.9 

1 204.5 

1 207.2 

1209.9 

i ■ 

0>u 

■J*0 

1 1 >A.O 

1 1 ^ U* •• 

1231.9 

1203. 7 

1 205.6 

1 208 . 1 

1210.8 

1213.5 

A 

v>* J 

V .0 

0*0 

1202.3 

1205.4 

1 207.5 

1 200. 3 

I Z1 1 .6 

1214.4 

1217.1 

5 

0*0 

G.O 

0.0 

0. 0 

1209.3 

1211.3 

1213.0 

121S.4 

1210.0 

1220.6 

t 


c • ^ 

J.O 

0. 0 

o.u 

12 t A. 1 

1215. 7 

1210.1 

1220. 7 

1223.3 


\ 1 

l£ 

13 

IH 

1 5 

1 b 

1 7 

10 

19 

20 

1 

1 >•<. 

1 « 1 1 • '« 

1 < 1 4*b 

121 7. 1 

1219.5 

1 22 1 .7 

1 223.2 

1 224.4 

1225.3 

1224.S 

2 

1 <: 1 ^ > 7 

1 £ 1 

1 £ 1 d • J 

1 22 J. b 

1 222 aU 

1 2 2 * . 1 

. 226 .5 

1227.7 

1 228.5 

1227.6 


l£ Ic .2 

1 e 12 .a 

1 1 *4 

1223. G 

Ic2t.2 

122c. 4 

1 229.9 

1231.0 

1 231. 9 

1231. 4 

4 

1 ’ t V • 7 

Mlimi 

t £24«<i 

1 ?2 7. J 

1 229•^ 

I2J 1 .8 

1 233.2 

1234.3 

1 235. J 

0.0 


J ■ 

1 _ ^ 

1 : r = « j 

1 j j : . 7 

l J .1 3 . 0 

1 ’ '5.2 

1236.6 

1 237.7 

0.0 

0.0 

V. 

•ilr / 

1 

i. 

1 C 1 -4 

1 c * «l 

J • u 
J . J 
■* 

i C £ C. • ^ 

1 

w » j 

J 

: ■ 7 

c 


I233i> 3 

1 2 a5.»> 

lt37.7 

1 239. 1 

0.)^ 

0. 0 

0. 0 


iNlEGKArno h£AT INPUT Ai.EACh NQO£. BIVJ 


HOW 

1 


PO« / CCL 


he* / cu. 


65.7 

201.7 

300.6 

327.6 

33 Ti % 

169.2 

ZT &.0 

266.4 

260.6 

256.2 

0.0 

2 C 4.6 

0.0 

0 . 0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3 J 6 .;£ 

0.0 

0.0 

0.0 

p.o 

0.0 

0*0 

0.0 

0.0 

0.0 

0^0 

353 . 6 

6.0 

0*0 

0.0 

0.0 

0.0 

0 * 0 . 

0.0 

0.0 

0 . 0 

0.0 

364*3 

0.0 

0*0 

0*0 

0 . 0 

0.0 

0*0 

0 . 0 ' 

0.0 

0.0 

0.0 

185.0 . 

307.3 

306.5 

306.4 

300.2 

1 1 

12 

J 3 

14 

15 

16 

1 7 

10 

19 

20 

250.2 

241 

23 C .2 

233.2 

227.9 

220,0 

85 . 4 

165.2 

155*0 

135.9 

0.0 

a.o 

0 . 0 

Q.O 

0,0 

0.0 

0*0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

' o.o 

0.0 

233.1 

J.O 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0 . 0 

0.0 

24 5.2 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

249.9 

0.0 

0.0 

JO 8 . t > 

307.3 

30 7 . 4 

309 . 0 

308,9 

306.6 

125.4 

C .0 

0.0 

0.0 


1 0 c* 5 
^1 1 • A 
J .0 
O.O 

a«o 

0 .0 


£2 

A »• J 

0.0 

0.0 

Q*0 

QmJ 

c.o 
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Section 5 


SPECIAL OUTPUT FILES 


5.1 DISCUSSION 

A option has been added to CAVE that will create a special file containing heat 
transfer coefficients, average adiabatic wall temperatures, and node temperature 
for pre-selected nodes at each print time. This file can be useful for making quick 
comparisons of these parameters in areas of the geometry that are critical. This 
file also can be used for generation of time-history plots of these parameters. The 
user can select this option for any of the available geometries. 

The user must assign a logical unit number for the special output file. The 
variable NODOUT, defined at line 1030 in Subroutine OVLY21, must be set equal to 
this logical unit number, and provision must be made in the user's job control 
language (JCL) to allocate a disk file, tape unit, etc. for this purpose. 

5. 2 INPUT DATA FORMAT 

The inputs that specify the nodes for which data will be written into the special 
file follow the Initial Temperature Cards. For example: 

Initial Temperature Cards 

• KODE, I, T(I), II, JJ 

• ... 

• ... 

• 11100 (signifies end of temperature cards) . 

Special Output File Card(s) 

• NPLOTS, ITCODE (214) 

• NODNUM(l) , NODNUM(2) , . . . NODNUM (NPLOTS) 

(Omit if NPLOTS = 0) (2014) 

NPLOTS = Number of nodes whose data will be written into 

special output file. If NPLOTS = 0, option is not 
being used. 
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ITCODE = 1 If node temperatures in special file are to be in R 

ITCODE = 2 If node temperatures in special file are to be in K 

NODNUM(I) = Node number of temperature to be stored 1=1, 2, 

3 . . . NPLOTS. 

Note : If the special output file option is not used, NPLOTS cannot be omitted 
and must be read in as zero. 

5. 3 SAMPLE INPUT & OUTPUT DATA 

The following pages show sample input data to exercise this option and the 
resultant output data. 
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SHEET 5.3 SAMPLE INPUT DATA FOR SPECIAL OUTPUT FILE (SHEET 1 OF 2) 



CJ1 


NODE NUMBERS 



(Q « 200^E i^ (0 ♦ 55^0E -02i (57i 50 WE:i04) CO > 5293. 4t.’i03) 
0.50000E+02 0»55500E --02 0 ♦ 15000E.- ♦ 04 0 *56230E' i03 


lOOOOE+03 

15000E+03 


•02 0*15000E-E04 
•02 0.15000E-f04 


0*61251EM03 

0*6'J694E-E03 


1827-028B 


0*20000E^f02 

0.0 

0*0 

0*51£J99E-f03 

0,50000E-f02 

0*0 

0*0 

0*55235EM-03 

0tl0000E-f03 

0.0 

0*0 

0*60274E-f03 

0»1E5000E^I03 

11 

0*20000E^I02 


0 * 

0*64738E-I03 

0»5^5500E-02 

o~^T5oo()F:Tt^ 


0»50000EM02 

Ot 55500 E--02 

0*15000E-f04 

o*55967;rt^ 

0»10000E-(03 

0 * 5 5500 E- 02 

0*15000E-(04 

0*A0947E-(03 

0.15000E+03 

0 » 55500E^^^^02 

0* 15000E-I04 

0*65362E-f03 

24 

0*20000E'f02 

0^0 

0*0 

0*51606E+03 

0»50000Ei02 

0*0 

0*0 

0 * 5486 7E •103 

o»iooooi;mo3 

OtO 

0*0 

0*59770E-I03 

0*15000E+03 

0*0 

0*0 

0*64137E-l-03 

29 

0*20000E-I02 

0»0 

0*0 

0*51469E^I03 

0*50000E^f02 

0*0 

0*0 

0*54705E-l-03 

0»10000E+03 

0*0 

0*0 

0*59567E-f03 

0*15000E-(03 

0*0 

0*0 

0*63907E’f03 

47 

0»20000E:M02 

0*0 

0*0 

0*51974E^»03 

0,50000E-f02 

0*0 

0*0 

0*54469EI03 

0»10000E-I03 

0*0 

0*0 

0*58426E-f03 

0»15000EI03 

0*0 

0*0 

0*62134E-i03 


NODE TEMPERATURE AT 
'this TIME 


.THESE VALUES ARE 
'ZERO BECAUSE THE 
NODES ARE NOT 
SURFACE NODES 
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Section 6 


CONVECTIVE HEAT TRANSFER COEFFICIENTS FOR LEADING-EDGE GEOMETRY 


6.1 DISCUSSION 

The CAVE code has been expanded to provide the user the following options for 
computing the heat transfer coefficients on the leading-edge geometry: 

• Option 1 - Entire flow over the leading-edge is assumed to be laminar 

• Option 2 - Entire flow over the leading-edge is assumed to be turbulent 

• Option 3 - Flow over the leading-edge is laminar up to a specified (i.e. , 
input value) location where the flow transitions is turbulent for the balance 
of the leading -edge. 

Options 1 and 2 were available in the original CAVE code. No change was made 
to these options. Option 3 is new and requires that the user specify the locations of 
the laminar to turbulent transitions on both the upper and lower leading -edge surfaces. 
Subsection 6.2 gives the required input data for the three options. 

For option 3, subroutine LEES2 calculates the turbulent heat transfer coefficient 
downstream of the transition point with suitable modification to reflect the upstream 
laminar flow influence on the boundary layer thickness. The development of the 
modification is now presented. 


The turbulent boundary layer thermal displacement thickness can be defined as: 
»*th = 


where C,p is a coefficient encompassing fluid properties and flow conditions, and S 
represents the distance measured along the surface from X = 0 to the point at which 
flow transition occurs. 


For laminar boundary layers, 3*^^ is proportional to the total heat transfer. 

^ _ ^laminar 

’ ® th , laminar V^ - H^) 


where H and H are the freestream and wall enthalpies. 
« w 
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Because Q. . = Qi 

^turb lam 


/ ^o, turb\ 
' ^o, lam / 


(h^ = stagnation point heat transfer 


coefficient) 

and because at the point where transition occurs the following relationship must be 
satisfied: 


^thjturb ^thjlam ’ 


it follows that 


I* = r c? 4/5 _ Sam 

"th.turb T eff 


Q 


turb 


p V (H -H ) 

00 00 00 W 


/^o,lam \ 
\^o,turb / 


(H^ - H^) 


where is the effective position of the transition point. 


However , 


C 

T transition 


Q . 
_ ^turb 


P^V„ (H^ - H ) 
00 00 ' 00 


therefore: 


/^o,lam \ 
' o,turb / 


^eff ^trans 


/h , v5/4 

/ o,lam \ 

X^Ojturb / 


This equation gives the relationship between the effective and the actual positions 
of the laminar to turbulent transition. 

/-s 

Within LEES2 the integral / P U dS (used to compute the heat transfer 

O"' 

coefficients) has been modified to account for transition. The following equation is 
employed for locations at or beyond the transition point: 

S .h , . 5/4 S, S 


r / 0 , 1 am \ ' i^trans C 

J j J PUdS J 

o ^ o.turb ' J S 


P U dS 
trans 


6.2 INPUT DATA FORMAT FOR LEADING-EDGE GEOMETRY WITH LAMINAR 
AND/OR TURBULENT HEATING 

The input data for the leading-edge geometry is unaltered with one exception. 
For the laminar to turbulent transition case, two additional input values are required: 
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STRA.NU and STRANL which specify the transition locations. The so-called Wing 
Angles card, part of the input data package for the leading-edge geometry, now 
includes these two new parameters. The format is as follows; 


Wing Angles Card 


SWEEPA, DIHEDA, CODEX, HMODI, TURBL, STRANU, STRANL (7E10.0) 

SWEEPA = wing sweep angle, in degrees 

DIHEDA = wing dihedral angle, in degrees 

CODEX = 0. for convective coefficient and adiabatic wall temperature 

computed by CAVE; = -1. for tabular input of coefficients 
and temperatures 

HMODI = nonzero value indicates that two tables will be read at 

the end and used to multiply the convective coefficient 

TURBL = less than zero, entire flow is assumed to be laminar 

= 0. entire flow is assumed to be turbulent 
= greater than zero, a location has been specified on the 

upper and lower leading-edge surface where flow transitions 
from laminar to turbulent 

STRANU = X location at which flow transitions on the upper surface 
(required only for TURBL >0) 

STRANL = X location at which flow transitions on the lower surface 

(required only for TURBL >0). 
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Section 7 


MODIFICATION TO RADIATION CALCULATION PROCEDURE 


Heat transfer problems involving radiation are non-linear and therefore not 
exactly tractable with the CAVE eigenvalue-eigenvector solution technique. The 
linearization procedure used in the original CAVE code is the common one of 
modifying the convection coupling to account for the radiation (see Appendix C of 
Reference 1). The temperatures of the surface nodes at the beginning of the time 
subinterval (end of previous time subinterval) established the radiation modifiers. 
Using temperatures at the beginning of the time subinterval is quite satisfactory for 
problems in which radiation plays a. minor role or where the surface temperatures 
are varying slowly from one time subinterval to the next (e.g. , "small" time 
subintervals) . 

The code has undergone modifications that permit larger time steps to be 
utilized in problems where radiation heat transfer is a significant factor. The 
modifications are such as to provide estimates of the surface temperatures at the 
end of the time subinterval which are then used with the temperatures at the beginning 
of the time subinterval to establish modified convection couplings for the time sub- 
interval (i.e. , the code now uses a mean temperature for the subinterval rather than 
the temperature at the start of the subinterval) . 

The following discussion outlines the new procedure that has been incorporated 
into the code. 

Consider a node diagram for a surface node of a body: 



1827-029B 
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where 


hA 

aeA 

kA 

L 

’’aw 



convection coupling, Btu/sec-R 
radiation coupling, Btu/sec-R^ 

conduction coupling, Btu/sec-R 

adiabatic wall temperature of fluid (more properly referred to as 
recovery temperature of fluid, R 

background radiation temperature (usually taken as 0 °R) , °R 
temperature of surface node, R 
temperature of interior node, R 


The heat transfer into the surface due to convection and radiation is 


This may be rewritten as: 


(Eq. 2) 


where 


Q,^=(h+hj,)A(T^^-T^) 

«(4-4) 

<^aw-V 


(Eq. 3) 


(Eq. 4) 


Eq. (3) can be rewritten into a form suitable for CAVE: 

^in ^ ’^eff ^ ^'^AW " V 

where 

hgff = h (E«l- 6) 

In the original version of CAVE, T^ was the temperature of the surface at 
time t, the start of the time interval. In the new version, T^ is the mean of the 
actual temperature at time t and the estimated temperature at t + At. This estimated 
temperature is obtained using the h^ values based on the T^ values for time t. 
Fundamentally, it is a five step process: 

1. Calculate h^^ based on T^(t) . 

2. Calculate T^ (t+ At) using eigenvalue-eigenvector procedure with hj^ 
from step 1, 
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3. Calculate mean of (t) and (t + At). 

4. Recalculate using mean value of from step 3. 

5. Recalculate temperature distribution in body using from step 4. 

It is noted that is typically negative (see Eq. (4) with > Tj^) and there- 
fore for conditions with the radiation coupling essentially equal to the convection 
coupling, h^jj approaches zero. If h^^^ approaches zero over the entire geometry, 
ill-conditioned matrices may result and the cumulative effect of arithmetic roundoff 
may preclude convergence of the eigenvector-eigenvalue iteration procedure within 
subroutine IJEN. A statement to the effect that convergence was not achieved will 
be printed by CAVE in such situations. The temperatures printed may however be 
fairly accurate in spite of the lack of convergence. This may be attributed to the 
eigenvector values being reasonable although not within tolerance and furthermore not 
being of great significance since the body is not changing very much in temperature 
due to the very small coupling between it and its environment. Double precision 
arithmetic would help to alleviate this "ill-conditioned" difficulty. 




Appendix 
FLOW CHARTS 


This appendix presents flow charts for the new subroutines and subroutines 
that have undergone more than trivial modification. 

All CAVE subroutines that have undergone any modification whatsoever have 
been given slight name changes to avoid any possiblity of the wrong subroutine being 
used. 


TABLE A-1 SUBROUTINES USED IN CAVE (SHEET 1 OF 2) 


Name 

Called By 

Main Purpose 

SIZES 

CAVE (main prog) 

Computes storage locations needed. Compares to 
number requested 

OVLYIO 

SIZES 

Sets up arrays for PCP5 

PCP5 

OVLYIO 

Reads & writes property data, controls geometry 

SLAB2 

PCP5 

Computes volumes & conduction shape factors for 
slab & cooled panel geometries 

X24C 

PCPS 

Computes volumes & conduction shape factors for 
basic X24C geometries 

GEN2 

PCP5 

Reads & writes volumes & conduction shape factors 
for general geometries 

MATOUT 

PCP5 

Writes material properties 

XTABS1 

GEN2 

Reads tabular values of hA & Taw Qeneral 

geometry problems 

OVLY21 

SIZES 

Reads initial temperature distribution & flight 
trajectory. Controls problem solution, steps time, 
computes average convection couplings for each step. 
Writes solution each time step 

UNFIT 

OVLY21 

Finds value from a table by linear interpolation 

ATTAC3 

OVLY21 

Finds node number closest to stagnation point & 
renumbers nodes as required by LEES2 

LEES2 

ATTACS 

Computes ratios h/h & Taw variation around 
leading edge problems 

NURED1 

OVLY21 

Reads tabular values of hA & TaW functions of 
distance & time 

DINTK 

OVLY21 

Finds value from a table using double interpolation 

PROP 

OVLY21 

Computes conduction douplings & mass specific heat 
product for each element given conduction shape 
factors & volumes 

PROUT2 

OVLY21 

Writes node numbers, conductances, capacitances 
& initial temperatures 

XINTPl 

OVLY21 

Finds values of several dependent variables from a 
table by linear interpolation on a single independent 
variable 

FLATH 

OVLY21 

Finds h & Taw over a flat plate 

ATMOS 

OVLY21 

Finds atmospheric pressure, temperature & density 
for given altitude 

POLRT 

FLATH 

Computes the roots of a polynominal 

TRANS 

FLATH 

Determines whether flow over flat plate will be 
considered laminar or turbulent 

DESDA2 

1827<030B 

OVLY21 

Calls eigenvalue & matrix routines. Calculates 
temperatures 



TABLE A-1 SUBROUTINES USED IN CAVE (SHEET 2 OF 2) 


Name 

Called By 

Main Purpose 

IJEN2 

DESDA2 

Obtaines dominant eigenvectors & eigenvalues of 
a given matrix (using Jennings method of simultaneous 
vector iteration) 

EIGVC 

DESDA2 

Prepares approximate guesses for the eigenvectors 
to start the Jennings algorithm iteration for the 
time step 

BFACS 

IJEN2 

Factorizes a banded positive-definite matrix 

BSOLS 

DESDA2, IJEN2 

Using the factors of a given banded positive-definite 
A as generated by BFACS solves for X in the system 
AX=Y 

ORNML 

IJEN2, DESDA2 

Carries out the standard Gram-Schmidt orthonormaliza- 
tlon of a group of vectors 

HETRAl 

DESDA2 

Sets up coefficient matrix (of conductances) In 
compact form 

RVORDR 

IJEN2 

Reorders estimated eigenvalues according to 
magnitude 

AORDER 

IJEN2, RVORDER 

Sets up permutation indices needed for ordering 
the eigenvalues 

DISPL2 

(Various) 

Prints Information, mainly debug special output. 
In array form 

PART 

DESDA2, PCP5 

Prints debug output information & intermediate 
timing of calculation 

DATE 

PART 

Determines data of run 

SWITCH 

DISPL2 

Converts columns of a matrix to rows of visa versa 

SCAPR2 

(Various) 

Computes scalar product of two vectors 

WALL2D 

SIZES 

Computes volumes & conduction shape factors for 
Scramjet wall geometry 

BLUNT2 

PCP5 

Computes volumes & conduction shape factors for 
solid & hollow (cooled) blunt nosed cone geometries 

CONSET 

OVLY21 

Sets up input for LEES2 for analytical solution of 
blunt nosed cones based on a flight trajectory table 

CHECKD 

(Various) 

Checks completion codes returned from Subroutine DINTK 

CHECKIM 

(Various) 

Checks completion codes returned from Subroutine NURED 

ROT 

POUT 

Prints headings for nodal array data for scramjet wall 
geometry 

POUT 

1827-031B 

WALL2D 

Prints out nodal array data for scramjet wall geometry 
In a format which resembles the geometric shape 
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SIZES: 


C 


START 


) 


i 


CAVE; ASSIGN TOTAL 
DIMENSION 



COMPUTE DIMENSIONS 
OF ARRAYS 




NURED1; READ IN STANDARD TABLE 
OF TEMPERATURE AS A FUNCTION OF 
ENTHALPY & LOG^q(P/Pq) 



1 


2 




LEADS: (LEADING-EDGE GEOMETRY) 

READ & WRITE MCAP, THETA, AX's, AY's, ETC 
COMPUTE CONDUCTION SHAPE FACTORS 
COMPUTE VOLUMES 

COMPUTE CONVECTION COUPLINGS FOR 
COOLED SURFACES 




X24C: (X24C STRUCTURAL ARRANGEMENT) 
READ & WRITE AX's, AY's, SI, S2, ETC 
COMPUTE CONDUCTION SHAPE FACTORS 
COMPUTE VOLUMES 

COMPUTE CONVECTION COUPLINGS FOR 
COOLED SURFACES 


40R 5 


BLUNT2: (BLUNT-NOSED CONE GEOMETRY) 
READ & WRITE AX's & AY's, ETC 
COMPUTE CONDUCTION SHAPE FACTORS 
COMPUTE VOLUMES 

COMPUTE CONVECTION COUPLINGS FOR 
COOLED SURFACES 


1827-033B 


FIG. A-1 ORGANIZATION OF CAVE IN TERMS OF THE MORE IMPORTANT 
SUBROUTINE CALLS (SHEET 1 OF 2) 
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FIG. A-1 ORGANIZATION OF CAVE IN TERMS OF THE MORE IMPORTANT 
SUBROUTINE CALLS (SHEET 2 OF 2) 
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C ENTER J 
INPUT: 

MCAP. THETA, 
AX, AR,TAU. 

RP. HCOOL, 
TCOOL 

~~r~ 


COMPUTE LENGTH 
OF CONE SECTION 
(AFT OF SPHERICAL 
NOSE) 


COMPUTE SPHERICAL 
NOSE RADIUS 



COMPUTE POSITION 
PARAMETERS FOR 
SPHERICAL NOSE 
SURFACE BOUNDARY 
NODES (PHI, EX, A) 


I 


COMPUTE POSITION 
PARAMETERS FOR 
AFT CONICAL 
SECTION 


I 


COMPUTER OUTER 
RADII OF SPHERICAL 
NOSE NODAL 
ELEMENTS 


I 


MOVE INNERMOST 
NODES TO CAVITY 
BOUNDARY & 
ADJUST VOLUMES 


COMPUTE VOLUMES 

nP ^PHFRIPAl 


Oi 1 

NOSE ELEMENTS 




COMPUTE VOLUMES OF 
INNERMOST NODES AT 
CENTER LINE 


COMPUTE VOLUMES 
OF RING ELEMENTS 
OF AFT CONICAL 
SECTION 



MOVE INNERMOST 
NODES TO CAVITY 
BOUNDARY & 
ADJUST VOLUMES 
OF RING ELEMENTS 


COMPUTE VOLUME OF 
INNERMOST NODES AT | 
CENTER LINE IN AFT 
CONICAL SECTION 


ASSIGN MATERIAL 
NUMBER TO ALL 
NODES 



1827-035B 

FIG. A-2 SUBROUTINE BLUNT2 FLOW CHART (SHEET 1 OF 2) 
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n 

i 
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FIG. A-2 SUBROUTINE BLUNT2 FLOW CHART (SHEET 2 OF 2) 


I 




ENTER 


/READ INPUT 
' DATA TO 
DEFINE SHAPa 
OF SCRAMJET/ 
WALL / 


PRINTOUT 
INPUT DATA 


CHECK ARRAY 
SIZES FOR 
THIS 

GEOMETRY 


SET ERROR 
FLAG 


/PRINTOUT , 
'array SIZEsi 
& ERROR 
MESSAGE / 


CHECK WALL 
DIMENSIONS 
FOR INPUT 
ERRORS 


RETURN 


COMPUTE AX'S 
FOR WALL'S 
FORWARD 
SECTION 


"dimensions 


SET ERROR 
FLAG 


PRINTOUT 
ERROR 
MESSAGE , 


COMPUTE NODE 

COMPUTE NODE 

COMPUTE A/X & 

VOLUMES FOR 



FORWARD 

AFT SECTION 

■ 

NODES 

SECTION 


i 


REORGANIZE 
AX'S IN CENTER 
SECTION FOR 
PROGRAM 
COMPATIBILITY 


COMPUTE NODE 
VOLUMES FOR 
CENTER SECTION 


COMPUTE NODE 
DISTANCE Z 
ALONG TOP OF 
WALL 


COMPUTE NODE 
DISTANCE Z 
ALONG BOTTOM 
EDGE OF WALL 


SET UP ARRAY 
IDENTIFYING 
BOUNDARY 
NODES (ID 
ARRAY) 


SET UP SURFACE 
NODE INDEXES 


COMPUTE AX'S 
FOR AFT 
SECTION 

1827-037B 


J /PRINTOUT / V 

/ AX'S, AY'S / 

/ & VOLUMES /' 

/ FOR ALL / 

/ NODES / 

FIG. A-3 SUBROUTINE WALL2D FLOW CHART 


RETURN 
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ENTER 



^ANTS NOD^ 
'N. TABLE 

I NO 


READ NODE 
NUMBERS FOR 
SPECIAL TABLE 


^GENERAL^ 

GEOMETRY. 

\ / 


y/^ SLAB X. 

/cooled panelTv 

X24C GEOMETRY, 

^cramjet wall, or 
\blunt cone^/ 


READ: 

REFX,C0DE X, 
HMODI 


SLAB,"V 

rcOOLED PANElS 
X^^X24C/^ 


LEADING-EDGE 

GEOMETRY 


& Taw^ 
IN TABLE 
s. FORM ^ 


SCRAMJET 
. WALL ^ 


& Tav^* 
IN TABLE 
s^^^FORM ^ 


BLUNT 

CONE 


^ &TAWX, 
IN TABLE 
sForm /• 

1 YES 


^CALL NURED 
JREAD TABLES) 
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FIG. A-4 SUBROUTINE OVLY21 FLOW CHART (SHEET 1 OF 6) 






FIG A-4 SUBROUTINE OVLY 21 FLOW CHART (SHEET 2 OF 6) 
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E 



1827-040B 

FIG. A-4 SUBROUTINE OVLY21 FLOW CHART (SHEET 3 OF 6) 


73 







H 


CALL LINFIT 
(GET TRAJECTORY 
PARAMETERS) 


CALL ATMOS 


CALL FLATH 
(COMPUTE h & 
Taw) 


COMPUTE AVG 

h & Taw 
FOR THIS 
TIME 


h & Taw > 

DISTRIBUTION 

■S^IFORM^ 




600 



105 

^/wERE^v 

COMPUTE 

INSERT h, hA, 

/ h&TAlA/ 

INDEXES FOR 

& Taw avg in 


TABLES? 

TOP & BOTTOM 

STACK ARRAY 

n 


SURFACE NODES 


MODIFY h FOR 

DIFFERENT 

LENGTHS 


COMPUTE 
INDEXES FOR 
BOTTOM 

SURFACE NODES 
ONLY 


^ IS \ 
FLOW 

TURBULENT 

V ? 


COMPUTE HEAT 
TRANSFER CO- 
EFFICIENT 
FOR TURBULENT 
FLOW 


CALL DINTK 
(INTERPOLATE 
h & Taw) 
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FIG. A-4 SUBROUTINE OVLY 21 FLOW CHART (SHEET 5 OF 6) 


75 














FIG. A-4 SUBROUTINE OVLY 21 FLOW CHART (SHEET 6 OF 6) 





77 


I 






1827-045B 


FIG. A-6 SUBROUTINE CONSET FLOW CHART 
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